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Abstract 


Three  kinds  of  rare  gas  (R)  containing  clusters  are  considered  in  this 
article:  (I)  ionic  clusters  R^  and  (R^M)*,  M  being  another  atom  or  a 

molecule,  (II)  neutral  rare  gas-halogen  clusters  RnX  with  R  -*•  X  electron 
transfer,  and  (III)  Rydberg  excited  clusters.  These  clusters  have  much  in 
common,  as  all  of  them  contain  ions  or  ionic  cores,  with  the  positive  charge 
delocalized  mostly  between  some  number  of  atoms.  Such  systems  can  be  treated 
by  the  semiempirical  diatomics-in-ionic-systems  (DIIS)  method.  We  present  In 
this  article  a  brief  description  of  the  DIIS  method  and  the  results  of  its 
applications  to  the  kinds  of  clusters  mentioned  above.  The  article  presents 
also  discussions  of  results  obtained  by  other  theoretical  methods  and 
experimental  studies. 
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I .  Introduction 


The  interaction  between  neutral  ground- state  rare  gas  atoms  is  of  a 
simple  nature  since  they  (i)  have  closed- shell  electronic  structure  and  are 
unable  to  form  chemical  bonds,  (ii)  are  spherically  symmetrical,  and  (Iii)  can 
be  excited  to  Rydberg  states  only.  Due  to  these  features,  the  structure  of 
the  molecules,  clusters,  liquids  and  solids  formed  by  rare  gas  atoms  is 
determined  mainly  by  the  atomic  sizes  and  the  interatomic  attractive  van  der 
Waals  (vdW)  forces.'*'’*’  The  interaction  becomes  more  complicated  when  other 
kinds  of  atoms  or  molecules  are  involved,  but  except  for  some  very  special 
cases,  this  interaction  does  not  significantly  effect  the  electronic  structure 
of  the  rare  gas  atoms.  The  features  of  a  rare  gas  atom  are  changed 

dramatically  if  an  electron  is  removed  from  a  closed  shell  as  the  result  of 
atomic  ionization  or  excitation.  The  rare  gas  atom  with  an  electron 
deficiency  in  the  outer  shell  is  chemically  active  and  is  then  able  to  form 
strongly-bound  molecules  and  clusters  whose  structure  is  more  intricate  than 
that  of  the  vdW  clusters.10’^  In  the  case  of  a  deficiency  on  an  inner  shell, 
the  atom  is  unstable,  which  leads  to  photon  emission  or  autoionization. 

We  can  distinguish  three  kinds  of  systems  formed  by  rare  gas  atoms  with 
electron  deficiency:  (i)  molecular  and  cluster  ions,  (ii)  molecules,  clusters 
and  solids  with  Rydberg-excited  rare  gas  atoms,  (iii)  molecules,  clusters  and 
solids  with  electron  transfer  from  a  rare  gas  atom  to  another  kind  of  atom-- 
halogen  atom,  mostly.  All  these  systems  have  much  in  common,  as  their 
properties  are  determined  to  a  great  degree  by  the  properties  of  the  rare  gas 
ionic  cores  which  are  the  same  in  the  ionic  and  excited  states.  The 
theoretical  consideration  of  these  systems  is  relatively  complex,  first  of  all 
because  of  the  presence  of  charged  particles  and  their  interaction  with  the 
neutral  atoms  which  they  have  polarized.  Further  complication  arises  because 
of  the  electrostatic  interaction  of  the  polarized  atoms  with  one  another. 
Dealing  with  the  systems  with  the  electron  deficiency,  we  also  have  to  take 
into  account  the  delocalization  of  the  electron  deficiency  between  rare  gas 
atoms.  Because  of  this  delocalization,  one  needs  to  take  a  quantum  chemistry 
approach  in  order  to  consider  such  systems.  Experimental  spectroscopic 
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studies  of  ionized  and  excited  molecules  and  clusters  containing  rare  gas 
atoms  have  much  in  common,  but  with  one  important  difference:  namely  that 
only  the  ionic  species  can  be  tested  by  mass  spectroscopy,  so  that  the  masses 
of  the  neutral,  both  ground  state  and  excited,  clusters  cannot  be  determined 
directly . ^ 

Considering  the  ionized  or  excited  clusters,  we  have  to  distinguish 
between  two  different  problems.  The  first  is  the  study  of  the  geometry  and 
spectroscopy  of  stable  ionic  or  metastable  excited  clusters.  Another  problem 
is  the  study  of  ionization  potentials  and  excitation  energies  of  neutral 
ground-state  clusters.  It  is  important  to  note  that  the  formation  of  ionic 
clusters  involves  often  as  an  intermediate  stage  the  Rydberg  excitations  whose 
study  is  necessary  to  understand  better  the  ionization  dynamics.  However,  the 
dynamical  problems  of  the  formation  and  decay  of  ionic  and  excited  clusters 
will  be  not  considered  at  length  in  this  review  article. 

Without  going  into  the  semantic  discussion  of  the  term  "cluster"  (there 
are  some  differences  on  this  point  in  the  literature),  we  shall  use  this  term 
here  to  refer  to  rare  gas  containing  systems  with  more  than  two  atoms  whose 
dissociation  or  atomic  detachment  energy  is  small,  less  than,  say  0.5  eV.  The 
name  "cluster"  shall  apply  not  only  to  stable  or  metastable  systems  but  also 
to  unstable  systems  formed  as  a  result  of  cluster  ionization  or  excitation. 

Any  stable  or  metastable  diatomic,  even  that  weakly  bound,  shall  be  called  by 
the  name  "molecule." 

In  the  next  section  we  shall  discuss  briefly  the  chemical  properties  of 
rare  gas  atoms  and  the  structure  of  neutral  vdW  clusters.  In  Section  III  we 
will  discuss  the  calculation  of  the  polarization  energy.  In  Section  IV  we 
will  review  the  semiempirical  method  called  diatomics - in- ionic  systems  (DIIS) 
which  we  have  developed  for  the  study  of  rare  gas  systems  with  charged  atoms. 
Ionic  clusters  will  be  considered  in  Sections  V  and  VI.  Excited  clusters  with 
electron  transfer  and  Rydberg-excited  clusters  will  be  considered  in  Sections 
VII  and  VIII,  respectively.  Concluding  remarks  are  given  in  Section  IX. 
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II .  Chemical  Properties  of  Neutral  Ground  State  Rare  Gas  Atoms 

In  their  neutral  ground  state,  '■are  gas  atoms  have  a  closed-shell 

C 

structure  with  the  (np)  outer  shell.  Because  of  such  structure  they  are 
chemically- inert  atoms,  at  least  as  they  interact  with  other  closed-shell 
atoms  or  molecules.  Being  chemically  inert,  rare  gas  atoms  interact  by  vdtf 
forces  only,  which  at  large  distances  are  represented  by  attractive  dispersion 
forces  roughly  proportional  to  the  polarizabilities  of  the  interacting 
species.  At  small  distances  the  dispersion  forces  are  overcome  by  strong 
repulsive  forces  which  determine  the  atomic  vdtf  radius.  The  rare  gas  atoms 
vdW  radii  and  polarizabilities,  as  well  as  the  ionization  potentials,  are 
presented  in  Table  1.  In  the  lowest  ionized  state,  the  outer  shell  structure 
becomes  (np)5,  with  an  electron  deficiency  on  a  np  orbital.  This  electronic 
structure,  similar  to  that  of  halogen  atoms,  makes  the  rare  gas  ions 
chemically  active.  The  diatomic  molecules  formed  by  rare  gas  ions,  such  as 
R+  10,18-20  and  21  - 27  ^  are  chemically  stable  molecules  whose 

dissociation  energy  exceeds  mostly  1  eV. 

The  geometrical  structures  of  the  weakly-bound  vdW  molecules  R,,  or  RR' 

and  clusters  R  are  determined  by  the  vdW  radii  of  their  component 
n 

atoiss.^'^'^'^  The  experimental  equilibrium  distances  and  the  dissociation 

29 

energies  of  the  rare  gas  diatomic  molecules  are  presented  in  Table  2.  The 

equilibrium  geometry  of  the  rare  gas  vdtf  clusters  can  be  described  as  the  most 

30 

compact  packing  of  almost  rigid  spheres.  For  example,  three  rare  gas  atoms 

form  an  equilateral  triangle  and  four  rare  gas  atoms  form  an  equilateral 

pyramid.  The  first  most  stable  Rn  cluster  with  magic  number  n-13  is  formed  by 

a  central  atom  and  12  atoms  symmetrically  located  around  the  central  one.  The 

30 

next  "magic"  cluster  R55  has  two  closed  coordinate  spheres  (1+12+42). 

Because  of  the  pairwise  character  of  the  vdtf  forces  and  their  short  radii,  the 

cluster  energy  is  practically  equal  to  the  sum  of  the  interaction  energies 

between  all  neighbor  pairs.  The  mutual  attraction  between  several  atoms  in 

vdtf  clusters  leads  to  some  small  decrease  in  the  equilibrium  interatomic 

distances  compared  to  those  of  diatomic  vdtf  molecules,  which  is  demonstrated 

for  Ar  clusters  in  Table  3,  where  the  dissociation  and  atom  detachment 
n 
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energies  of  these  clusters  are  presented  as  well.  It  is  important,  however, 

to  note  that  due  to  thermal  motion  in  the  weakly-bound  vdW  clusters,  the  atoms 

are  mostly  localized  far  away  from  the  equilibrium  (minimum  energy)  positions. 

P.are  gas  atoms  form  typical  vdW  clusters  also  with  nonpolar  molecules,  such  as 

31  33  35 

halogen  Xj  molecules'"  "  or  aromatic  molecules.  *  Somewhat  more 

complicated  and  more  strongly-bound  clusters  are  formed  by  rare  gas  atoms  with 

polar  molecules  which  contribute  to  attractive  polarization  forces  in  addition 

36-39 

to  usual  vdW  forces.  *  All  these  clusters  have  much  in  common,  as  the 
electronic  structures  of  the  component  atoms  and  molecules  are  unaffected  by 
one  another  (except  for  polarization  in  the  presence  of  polar  molecules) .  The 
mutual  independence  of  atomic  and  molecular  electronic  structures  in  vdW 
clusters  is  well  demonstrated  by  the  weak  effect  of  the  constituent  rare  gas 
atoms  on  the  valence  spectrum  of  a  molecular  cluster . Since  the 
electrons  in  the  vdW  clusters  are  strongly  localized,  the  cluster  wave 
function  can  be  presented  as  the  antisymmetrized  product  of  atomic  and 
molecular  wave  functions 


*  “  A  n  Xa 


(II. 1) 


where  Xj  are  the  group  functions  describing  individual  atoms  and  molecules. 

The  electronic  structure  of  rare  gas  atoms  can  be  affected  by  open-shell 
radicals  with  high  electron  affinity,  particularly  by  halogen  atoms  X.  Due  to 
the  high  electron  affinity  of  halogen  atoms,  the  excitation  energy  of  the 
states  with  electron  transfer  (R+X  )  is  relatively  low,  for  example  4  eV  in 
XeCl  molecule. ^  Taking  into  account  the  coupling  between  the  nonpolar  (RX) 
and  polar  (R+X)  electronic  structures,  the  RX  wave  function  has  to  be 
described  as  the  supperposition 


9  -  C1$fLX  +  C2$R+x'  (II. 2) 

Because  of  the  coupling  between  these  two  electronic  structures,  the  ground- 
state  R^X  systems  cannot  be  considered  as  pure  vdW  complexes,  since  the 
coupling  contributes  to  some  valence  attraction  between  the  halogen  and  rare 
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gas  atoms. ^  For  example,  the  contribution  of  the  coupling  (valence)  energy 
constitute  45X  of  the  XeCl  dissociation  energy,  and  the  Xe-Cl  equilibrium 

o 

distance  (3.23A)  is  much  smaller  than  the  sum  of  the  vdW  radii  of  Xe  and  Cl 

O  Q 

atoms  (4A) .  The  coupling  between  the  nonpolar  and  polar  (ionic) 
configurations  becomes  much  more  important  in  the  case  of  the  heavy  rare  gas 
atoms  Xe  and  Kr  and  light  halogen  atom  F,  leading,  e.g.,  to  the  formation  of 
XeF^  and  XeF^  molecules . 


III .  Polarization  Energy  and  Charge  Motion  Effects 

The  presence  of  charged  particles,  such  as  R+  ions  in  ionic  R*  clusters, 

ionic  cores  and  electrons  in  Rydberg  excited  clusters,  and  R+  and  X*  ions  in 

excited  clusters  with  electron  transfer,  results  in  the  polarization  of 

neutral  atoms  and  molecules.  When  the  charges  in  a  cluster  are  strongly 

localized  on  some  atoms,  the  polarization  energy  can  be  easily  represented  in 

49 

the  point-dipole  approximation  by 

p  -  p<D  +  p<2)  (III. 1) 


,(D  - 


•iAj  a>izij 

1  J  L  J 


(III. 2) 


>(2)  _  2 


e  X  aj  X,ay  E  rjj'  ,mn  £  qiZij,mZij'  ,n  ’ 
j  J  j  J  mn  J  J  l  J  J 


j'>j 


(III. 3) 


where  are  the  ion  charges  in  the  electron  e  units,  Oj  are  the  neutral  atoms 
polarizabilities  (Table  1),  and  A^  are  the  parameters  which  determine  the 
atomic  polarization  energy, 


Aj  -  °je  n 


(III. 4) 


.  are  the  vectors  proportional  to  the  electrostatic  field, 


Zij  -  (Ri  - 


(III. 5) 
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and  Rj  being  che  radius  vectors  of  charged  and  neutral  atoms,  respectively, 

and  T..,  are  the  tensors  of  dipole-dipole  interaction, 
jj  ,mn 


rjy„n-  3<Rj  ■  WRj  • 


(III. 6) 


where  m,n  -  1,2,3  denote  the  x,y,z  projections  and  <5^  is  the  Kronecker  delta. 
The  number  of  charged  particles  (ions)  is  I  -  1  in  ionic  clusters  and  1-2 
(X*  and  R+)  in  clusters  with  electron  transfer  from  R  to  X.  The  first  term, 

in  Eq.  (III.l)  gives  the  energy  of  the  interaction  between  the  charged 
particles  (i  -  1...I)  and  the  neutral  atoms  which  they  polarize,  whereas  the 
second  term,  ,  gives  the  energy  of  che  dipole-dipole  interaction  between 

the  polarized  atoms.  The  contribution  of  the  dipole-dipole  interaction  is 
supposed  to  be  small  compared  with  the  direct  charge-neutral  atoms 
interaction. 

The  polarization  energy  presents  some  difficulties  in  the  theoretical 
consideration  of  clusters  with  charged  atoms.  First,  the  polarization  energy 
is  not  of  pairwise  character,  except  for  the  case  of  one  charge  and  the 
neglect  of  the  second  (dipole-dipole)  term  in  Eq.  (III.l).  Second,  using  the 
electrostatic  expressions  in  the  point-dipole  approximation,  we  are  neglecting 
the  overlap  of  the  ion  and  the  atom  which  it  polarizes,  leading  to  the 
unphysical  divergence  of  the  interaction  energy  at  small  separations.  Third, 
the  ionic  charge  is  often  delocalized  between  two  or  more  atoms.  The  first  of 
these  problems  will  be  considered  in  the  next  section  within  the  framework  of 
the  DI IS  method.  The  unphysical  behavior  of  the  ion-atom  polarization  energy 
(III. 2)  at  small  separations  is  usually  avoided  by  introducing  a  damping 
function  7  into  the  polarization  energy  expression,5® 


PiA 


(III. 7) 


where  R  is  the  separation  between  the  ion  i  and  atom  A.  Different  damping 
function  expressions  have  been  proposed  for  the  ion-atom  ’  and  electron- 

atom54*56  interactions.  If  the  correct  asymptotic  conditions  are  satisfied, 
the  form  of  the  dependence  of  7  on  R  does  not  imply  much  about  the  results. 
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For  example,  according  to  our  estimates,  different  reasonable  models  for  the 

damping  function  7(R)  result  in  variations  of  the  energy  of  Ar4  clusters 

14 

within  a  range  of  only  0.003  eV.  We  have  suggested  to  take  the  damping 
function  in  the  form 


'(R)  -  1  +  C 


p.+p,  12 


(III. 8) 


where  p.  and  p,  are  the  atomic  radii  of  the  ion  and  atom,  respectively.  For 
L  A 

the  values  p.  and  p.  in  the  Eq .  (III. 8),  we  are  using  the  atomic  vdW  radii. 
When  the  interatomic  separation  R  is  larger  than  the  sum  of  atomic  radii  (non¬ 
overlapping  atoms),  the  function  7  is  close  to  unity,  i.e.,  the  point-dipole 
approximation  is  valid.  When  the  interatomic  separation  is  smaller  than  the 
sum  of  atomic  radii  (overlapping  atoms),  the  damping  function  decreases  as 
R'4,  so  that  the  ion-atom  polarization  energy  (III. 7)  is  saturated  with  a 
limiting  value  of  -(^/(p^+p^^ .  Introducing  the  damping  function  (III. 8)  into 
Eqs.  (III.  2-3),  we  obtain  the  following  expression  for  the  polarization 
energy  of  a  polyatomic  system: 


,(1)  -  -l  Ai  <£  *i  zij)2’  T'ii  - 


(III. 9) 


>(2)  -  e2  7  a.  V  a . ,  Y  T .  . ,  Y  q.  /y77  Z. . 

t  V  J  Sn  JJ  'mn  i  1  ^ 

j  }  mn  1 


In  the  case  of  one  charge  (I  -  1),  Eqs.  (III. 9-10)  become 


p(D  _  V  A  -y  z2 
P  »  - L  Aj^ijzij 


(III. 9' ) 


P«>  la  l  r 

i  j  j  mn  J  J 
J  j'Oj)  ran 


/y. .  Z..  /y . . ,  Z .  . , 
mn  lj ,m  rij  ij  ,  n 


(III. 10' ) 


The  description  of  the  charge  delocalization  between  two  or  more  atoms  is 
ilicated  bv  the  dynamics  of  the  charge  motion.^2  Let  us  consider,  for 
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example,  a  rare  gas  cluster  ion  R+.  In  small  R*  clusters  the  charge  is 
delocalized  between  three  atoms  which  constitute  the  ionic  core  of  the  cluster 
(for  a  discussion  of  the  structure  of  the  Rn  clusters,  see  Section  V),  whereas 
other  R  atoms  are  left  netural  or  almost  neutral.  In  many  other  ionic 
clusters,  such  as  heteronuclear  rare  gas  clusters  Ar2Nen  or  aromatic  cations 
with  rare  gas  atoms  M+Rn,^  it  is  also  possible  to  distinguish  between  the 
ionic  core  and  neutral  atoms  which  are  polarized  by  the  ionic  core.  The 
charge  delocalization  in  a  positive  ionic  cluster  core  (molecule)  can  be 
described  as  a  motion  of  a  positive  hole  with  some  characteristic  time  of  the 
hole  "jump"  between  adjacent  atoms.  This  time  is  on  the  order  of  h/E,  where  E 
is  the  coupling  energy  responsible  for  charge  delocalization.  Since  neutral 
atoms  and/or  molecules  also  need  some  time  to  be  polarized,  the  interaction 
between  them  and  the  delocalized  charge  is  influenced  by  dynamical  effects. 
When  the  polarization  time  is  much  smaller  than  the  characteristic  time  of  the 
hole  motion  (adiabatic  polarization),  then  the  polarization  will  follow  the 
hole  motion  and  the  polarization  energy  will  become  an  average  of  the 
polarization  energies  of  all  hole  localizations.  Thus  in  the  case  of  one 
neutral  atom  polarized  by  an  ionic  molecule  or  cluster  core,  the  adiabatic 
polarization  in  the  point-dipole  approximatation  is 

Pa  “  \  at  l  a?F?  ,  (III. 11) 

2 

where  aj.  is  the  hole  population  at  the  i-th  atom  and  is  the  field  generated 
by  a  hole  located  at  the  i-th  atom.  In  the  opposite  case,  where  the 
polarization  time  is  much  larger  than  the  characteristic  time  of  hole  motion, 
(diabatic  polarization),  the  neutral  atom  will  be  polarized  by  the  average 
field, 


(III. 12) 


As  an  example  of  a  system  (unrealistic,  unfortunately)  which  demonstrates  a 
dramatic  difference  between  two  kinds  of  polarization,  we  suggest  an  ionic 
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ring  with  a  neutral  atom  at  the  ring  center.  For  the  diabatic  case  it  is 

4 

zero,  whereas  for  the  adiabatic  case  it  is  -C7(R)/R  ,  R  being  the  ring  radius. 

In  the  simple  case  of  a  diatomic  ion  and  the  neutral  atom  which  it 
polarizes,  the  general  (from  a  dynamical  point  of  view)  expression  for  the 
polarization  energy  is 


P  -  - 


4(e  +  2h)  f£(7iFi  +  '12F2) 


h(i/  72Fi  +  v'7oF 


2^2)2J 


(III. 13) 


where  €  is  the  atomic  S  -»  P  excitation  energy  which  determines  the  atom 
polarization  time,  and  h  is  the  interatomic  exchange  energy  of  the  hole  motion 
which  determines  the  hole  motion  characteristic  time.  According  to 
Eq.  (III. 13)  the  adiabatic  (e  >  >  h)  and  the  diabatic  (e  <  <  h)  polarization 
energies  are 

pa  -  -  5  a  (7lfJ  +  72F%)  (III. 14) 

pd  -  -  5  a  +  ^F^)2  ■  (III. 15) 


The  transition  from  the  adiabatic  to  diabatic  polarization  takes  place  at 
h  =  0.5e,  so  that  the  more  accurate  polarization  criteria  are 


n  >  >  2h.  adiabatic 


(III. 16) 


«  <  <  2h,  diabatic  (III. 16' ) 

It  is  reasonable  to  assume  that  in  more  complicated  systems  these  criteria 
will  be  valid  as  well,  at  least  roughly. 

Dealing  with  ionic  clusters,  we  need  to  know  whether  to  use  the  adiabatic 
(III.  11)  or  diabatic  (III. 12)  expressions.  For  most  ionic  clusters,  the 
adiabatic  condition  (III. 16)  is  fulfilled,  since  the  exchange  term  h  is 
usuallv  on  the  order  of  1-2  eV ,  whereas  the  excitation  energy  e  is  generally 
on  the  order  of  10  eV.  In  the  literature,  however,  the  model  calculations  are 
usuallv  performed  within  the  diabatic  approximation  i.e.,  considering  the 
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51  53  58 

polarized  atom  in  the  average  field  of  a  delocalized  hole.  ’  ’  According 

to  our  estimations  for  Xe^R  (R*He,Ne,Ar)  and  benzene+-R  cluster.1:  the 
difference  between  the  adiabatic  and  diabatic  polarization  energies  is 
noticeable,  but  not  large,  accounting  for  some  10X  of  the  polarization  energy. 
The  adiabatic  polarization  is  expected  to  take  place  also  in  the  excited 

systems  with  electron  transfer,  such  as  the  unstable  Cl  Xe+Xen  complexes  which 

5  9 

are  formed  in  irradiated  Xe  solids.  In  these  complexes  the  pcsitively- 
charged  hole  is  delocalized  between  some  number  of  Xe  atoms,  like  in  the  ionic 
Xe^  clusters.  ^  The  polarization  of  rare  gas  atoms  by  tha  ionic  cores  of  the 
Rydberg  excited  clusters  has  to  be  considered  in  the  adiabatic  approach  as 
well.  However,  as  will  be  shown  in  Section  V,  the  polarization  of  the  rare 
gas  atoms  by  the  Rydberg  excited  electron  should  be  considered  in  the  diabatic 
approach . 

IV .  Diatomics-In-Ionic-Svstems  (DUS')  Method 

As  mentioned  in  the  previous  section,  the  ionic  charge  (hole)  in  the 
clusters  with  charged  rare  gas  atoms  is  mostly  delocalized  between  two  or  more 
atoms.  In  order  to  study  such  systems  theoretically,  we  have  previously 

g 

developed  the  semiempirical  diatomics- in- ionic-systems  (DIIS)  method.  This 
method  makes  use  of  diatomic  potentials  as  initial  data  for  a  calculation,  but 
in  a  different  way  than  the  diatomics  - in-molecules  (DIM)  method. 

Let  us  consider  a  polyatomic  system  (A^^-.-Ay)*  consisting  of  J  closed- 
shell  atoms  or  ions  A^  with  one  electron  deficiency  (hole)  delocalized, 
generally  speaking,  between  these  atoms  (ions).  For  example,  the  ionic  R* 
cluster  is  considered  as  (RR. . ,R)  +  and  the  neutral  (ground-state  and  excited, 
with  electron  transfer)  rare  gas-halogen  cluster  RnX  is  considered  as 
(RR...RX  )+.  The  spin  of  such  systems  with  one  electron  deficiency  is  S  -  j. 
The  DIIS  wave  function  is  presented  as  a  linear  combination  of  diabatic 
polyatomic  wave  functions  with  a  localized  hole  as 

J 

♦  -  l  l  ■ 

i-1  m-1 


(IV. 1) 
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where  $>.  Is  the  diabatic  wave  function  of  the  configuration  A,A~...A;  ...A,, 
with  the  index  m  indicating  the  orientation  of  the  At  open  shell  (M^  -  3  for  a 
P  symmetry  shell).  For  example,  the  rare  gas  ionic  cluster  R*  is  described 
by  3n  diabacic  configurati  ns  R...R+...R  for  m  -  1,2,3.  The  R£X  system  is 
described  by  nine  diabatic  configurations:  R+RX  ,  RR+X  and  RRX  (m  -  r,2,3). 

An  ionic  cluster  consisting  of  n  rare  gas  atoms  and  one  hydrogen  atom 
(RnH)+  can  also  be  treated  by  the  DIIS  method  if  the  electron  affinity  of  H  is 
neglected  and  singlet  (S  -  0)  states  only  are  considered.  The  (R^H)*  cluster 
is  described  by  3n+l  diabatic  configurations:  R...R  . . .RH  (M^  -  3)  and 
R. . .R. . . RH+  (M±  -  1) . 

The  diabatic  polyatomic  wave  function  is  presented  in  the  DIIS  method 

^  3 

as  an  antisvmmetrized  product  of  atomic  group  functions,  which  are  assumed 
to  be  exact  many-electronic  wave  functions  of  individual  atoms, 

-  J 

*im  ’  A  (J*j  W  •  J  **  1  <IV-2> 

A 

where  A  is  the  operator  which  antisymmetrizes  electrons  belonging  to  different 

group  functions,  like  in  Eq.(II.l).  Since  the  nonorthogonality  of  the  atomic 

group  functions  complicates  significantly  the  energy  expression, ®  the  zero 

overlap  of  atomic  orbitals  (ZOAO)  approximation  is  usually  used  in 

62  63 

semiempirical  methods,  in  particular  in  the  DIM  method.  ’  The  ZOAO 

A 

approximation  allows  one  to  omit  the  antisymmetrization  operator  A,  thus 
rendering  the  polyatomic  wave  function  (IV. 2)  as  a  simple  product  of  atomic 
group  functions.  Due  to  the  lack  of  interatomic  electron  permutation  in  the 
polyatomic  function  (IV. 2)  the  Hamiltonian  of  the  system  can  be  partitioned 
into  atomic  and  diatomic  terms . ® 

Th«  diagonal  matrix  elements  of  the  DIIS  wave  function  (IV. 1)  have  the 
physical  interpretation  as  the  energies  of  the  diabatic  configurations,  i.e., 
configurations  with  fixed  localization  of  the  hole.  Assuming  the  adiabatic 
polarization  to  be  valid  (see  Section  III),  we  include  the  polarization  energy 
(111.1,9-10)  into  the  diagonal  matrix  elements.  In  the  ionic  clusters  with 
one  charge  and  negligible  small  dipo le - dipole  interaction  (111.10'),  the 
polarization  energy  (III. 9')  is  expressed  in  a  pairwise  way  and  can  be 
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included  in  che  diatomic  ion-atom  interactions.  If  all  interactions  are  of 
pairwise  character,  then  che  diagonal  matrix  elements  are  presented  as  a  sun 
of  atomic  energies  and  diatomic  interactions  as 


U 

im,  im 


*  i 


(IV. 3) 


where  W.  is  the  energy  of  the  closed- shell  atom  A,  ,  W.  .  is  the  diatomic 
J  J  J  lJ  2 

potential  between  two  closed-shell  atoms  A.  -  A.  ,  U.  is  the  energy  of  the 

J 1  J  2  1 

open-shell  atom  At  with  an  electron  deficiency,  and  U^m  ^  is  the  diabatic 
diatomic  potential  between  the  m-oriented  open-shell  atom  A^m  and  the  ciosed- 
sheLl  atom  A y  In  the  case  of  ionic  clusters,  the  closed-shell  atoms  Aj  are 
neutral  and  the  open- shell  atoms  At  are  the  ions.  However,  in  other  systems 
the  open- shell  atom  with  electron  deficiency  At  might  not  be  an  ion  but 


rather  a  neutral  atom,  such  as  (Cl  )+  *  Cl 


In  clusters  with  more  than  one  charge,  such  as  those  with  R+  and  X'  ions, 
Rydberg-excited  clusters  (positive  atomic  core  and  excited  electron)  or 
multiply-charged  ionic  clusters , as  well  as  in  any  systems  with  non- 
negligible  dipole-dipole  interactions,  the  polarization  energy  cannot  be 
expressed  in  a  pairwise  way.  In  this  case  the  polarization  energy  has  to  be 
singled  out  into  a  separate  term  P^ ,  and  the  DIIS  diagonal  matrix  elements  are 
expressed  as  a  sum  of  and  diatomic  diabatic  (Coulombic)  potentials  0  as 


h.  .  *  y  w.  +  y  y  u.  .  +  u,  +  y  u .  .  +  p . ,  j » J  i » j  o  **  i 

Mim,  im  4*  j  j2  (>jx)  ^2  1  j  lm’J  1  Jl  2 


(IV. 4) 


where  0-  ■  are  the  Coulombic  potentials  between  pairs  of  closed  shell  atoms, 

J  1J  2 

and  0.  ■  are  the  Coulombic  potentials  between  open-shell  atoms  (with  a  hole) 

im,  j 

and  closed-shell  atoms.  The  ion-atom  Coulombic  potentials  are  obtained  by 

excluding  the  polarization  term  (III.  7)  from  the  potentials  W.  .  and  diabatic 

J1J2 


potentials  U-  • , 
“  im ,  j 


0.  i  -  w.  -  A.  7(R,  i  )/R,  i 

J 1-1 2  J1J2  J2  J1J2  J1J2 


(IV. 5) 
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0.  .  -  u. 

im,j  im,j 


V(RiJ)/RT».j 


(IV. 6) 


In  Eq.  (IV. 5)  the  j^-th  atom  is  an  ion,  and  in  Eq.  (IV. 6)  the  i-th  atom  is  an 

ion.  If  both  atoms  are  neutral,  the  Coulombic  potentials  coincide  with  the 

diatomic  potentials  W  and  U.  The  polarization  term  is  determined  by 

Eqs.  (III. 9-10)  with  charges  q^  -  1  (the  index  i  in  Eqs.  (III. 9-10)  is  for 

charged  atoms,  whereas  in  this  section  it  is  for  the  hole  localization). 

The  off-diagonal  matrix  elements  H.  -  between  the  diabatic  states 
°  im^ , lmj 

with  the  same  hole  localization  i  but  different  orientation  m  are  expressed  by 
the  diabatic  diatomic  potentials  j ,  whereas  the  off-diagonal  elements 
between  the  diabatic  states  with  different  hole  localizations  are  expressed  by 

g 

the  diatomic  exchange  terms  V 


J .  .  -IdVY.Y.  h..  V.  Y. 

1lml’i2m2  J  x2m2  1112 


(IV. 7) 


where  dv  is  a  diffemtial  volume  element,  h.  ■  is  the  Hamiltonian  of  the 

1112 

diatomic  interaction,  and  x ^  and  Xim  are  the  group  functions  of  the  closed- 

shell  and  open-shell  atoms,  respectively.  The  transformation  of  the 

potentials  U  and  exchange  terms  V  between  m- oriented  orbitals  to  the  2  and  II 

potentials  and  exchange  terms  is  described  in  Reference  8. 

The  diatomic  potentials  W.  .  between  closed-shell  atoms  are  taken  from 

J  lJ  2 

empirical  studies  or  ab  initio  calculations.  The  diabatic  potentials  U^m  j 

between  open-shell  (i)  and  closed  shell  ( j )  atoms,  such  as  R+-R,  have  to  be 

calculated  since  only  adiabatic  potentials,  such  as  W  are  the  physically 

R2 

defined  values.  The  calculation  of  the  diabatic  potentials  can  be  performed 
by  considering  the  diatomic  fragments  within  the  DI1S  approximation  and 
solving  the  inverse  problem  of  the  2x2  matrix  eigenvalues,  i.e.,  calculating 
the  matrix  elements  for  given  eigenvalues.  In  the  case  of  the  homonuclear 
fragment  R^,  the  2x2  matrix  of  the  fragment  is 

U-E  V 

-  n  (iv. 8) 

V  U-E 
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where  U  is  the  R+  R  diabetic  potential  and  V  is  the  exchange  term  in  the 

(R+-R) - (R-R+)  coupling.  Substituting  into  Eq.  (IV. 8)  the  known  adiabatic 

potentials  -  ^2  and  E^  -  ^Z  ,  we  can  easily  find  the  diabatic 

n  g 

potential  U  and  the  exchange  term  V  for  2  symmetry.  In  the  same  way,  U  and  V 
of  n  symmetry  are  found.  For  the  case  of  heteroatomic  fragments  there  are  two 
different  diabatic  potentials,  so  that  one  needs  an  extra  empirical  or 
ab  initio  value  to  solve  the  inverse  eigenvalue  problem  and  to  find  the 
diabatic  potentials  U  and  the  exchange  term  V.  For  such  an  extra  value  it  is 
possible  to  use  the  diatomic  static  or  transition  dipole  moment.  ’  After 
the  diabatic  potentials  U  and  the  exchange  terms  V  of  all  diatomics  fragments 
are  determined,  the  Hamiltonian  matrix  elements  (presented  in  Reference  8)  can 
be  calculated,  and  the  wave  function  (IV. 1)  and  its  energy  spectrum  can  be 
found. 

In  the  wave  function  (IV. 1),  the  spin-orbit  coupling  is  not  taken  into 
account,  at  least  directly.  However,  the  energy  level  shifts  resulting  from 
this  coupling  can  be  incorporated  into  the  system  energy  indirectly  by  taking 
as  the  input  of  the  DIIS  calculation  the  empirical  or  ab  initio  potentials 
with  spin-orbit  coupling.  By  presenting  the  wave  function  (IV. 1)  without 
spin-orbit  coupling,  we  are  losing  some  of  the  excited  states  and  are 
excluding  from  our  consideration  the  S-P  mixing  effect  in  the  collinear 
configuration. 

V.  Ionic  Rare  Gas  Clusters 

The  ionic  rare  gas  clusters  R*  have  been  studied  extensively  both 
experimentally  and  theoretically,  particularly  recently.  In  spite  of  this, 
the  structure  of  these  clusters  is  still  not  fully  clear.  Even  for  the  case 
of  the  smallest  clusters,  R^,  the  experimental  evidence  does  not  allow  one  to 
chose  definitely  the  electronic  and  geometric  structure.  From  a  theoretical 
point  of  view  there  are  two  alternative  structures  for  R^,  namely  an 
asymmetrical  RjR  structure  with  the  charge  concentrated  on  Rj  and  a  structure 
with  the  charge  distributed  among  all  three  atoms.  For  the  case  of  R^R,  where 
the  neutral  atom  is  bound  to  the  Rj  molecule  by  polarization  forces 
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predominately ,  the  triangle  geometry  is  expected  to  be  the  most  stable  one. 
This  triangle  geometry  was  supported  by  the  first  quantum  chemistry 
calculation  of  Ar^  performed  by  the  approximate  XQ  method. ^  When  the  charge 
is  delocalized  among  all  three  atoms,  the  exchange  interaction  contributes 
significantly  to  the  binding,  causing  the  symmetrical  linear  geometry  to  be 
the  most  stable  one.  Such  geometry  is  supported  by  relatively  accurate  Cl 
ab  initio  calculations  of  He^ ,  Ne2  and  Ar^ ,  as  well  as  by 

DIM70'71  and  other  semiempirical^  ’  7^  ’  7^  calculations  of  Ar^  ^’70-72  an(^ 

Xe^ . 7^ ’ 7^  In  contrast  to  these  findings,  a  recent  Ar^  ab  initio 
calculation  shows  an  asymmetrical  linear  Ar^Ar  geometry  to  be  the  most  stable 
one.74  Such  geometry,  although  a  little  strange  from  the  point  of  view  of  an 
energy  balance,  can  be  explained  by  some  contribution  of  the  exchange 
interaction  between  the  Ar^  dimer  and  the  slightly  charged  (-  0.1)  outer  Ar 
atom. 

The  experimental  data,  obtained  mainly  for  the  Ar^  cluster,  cannot 

provide  direct  evidence  about  the  geometry.  The  conclusions  made  by  the 

indirect  analysis  of  the  experimental,  mostly  spectroscopic,  data  are 

contradictory .  Whereas  most  of  the  experimental  work  confirms  the  symmetrical 

linear  Ar^  structure ,  ^  •  75-81  oC^er  substantiate  the  asymmetrical  Ar^Ar 

82  83 

structure,  either  of  triangular  or  linear  geometry.  Such  discrepancy  is 

due,  most  probably,  to  the  features  of  the  Ar^  potential  energy  surface 

corresponding  to  a  very  flat  Ar^  potential  in  the  collinear 
72  74 

configuration.  ’  Such  a  potential  causes  the  vibrational  motion  to  shift, 
for  the  most  part,  the  Ar^  configuration  far  away  from  the  equilibrium 

84 

configuration,  thus  significantly  affecting  the  electronic  transitions. 

This  problem  has  been  considered  in  detail  in  Reference  72  by  means  of  a 
trajectory  study.  According  to  this  study  the  symmetrical  linear  Ar^  minimum 
energy  configuration  alone  provides  a  spectrum  close  to  the  experimental  one, 
if  one  takes  into  account  the  vibrational  motion.  The  experimental  studies  of 
other  R^  clusters,  such  as  He2,  Kr^,  and  Xe^,  ’  were  concentrated 
mainly  on  the  problems  of  cluster  stability  and  dynamics  of  formation,  without 
much  consideration  of  the  structure.  The  interpretation  of  the  Xe^  spectrum 
obtained  in  Reference  88  supports  the  linear  geometry  of  this  cluster. 
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In  our  calculations^  performed  by  the  semiempirical  DIIS  method  (see 

4.  -f 

Section  IV),  the  triatomic  clusters  Ar^  and  Xe2  have  a  symmetrical  linear 
equilibrium  configuration,  in  accordance  with  most  of  the  other 
calculations . ^ ^  For  the  interatomic  distances  at  the  Ar^  equilibrium 

o 

configuration,  our  calculation  gives  2.59  A,  exactly  the  same  as  another 

72 

semiempirical  calculation  but  slightly  less  than  an  ab  initio  calculation 
°  69 

(2.62  A).  The  cluster  charge  is  shared  equally  between  the  inner  and  two 

2  + 

outer  atoms  (Figure  1),  like  in  Reference  72.  In  the  excited  II  state  the 

§ 

charge  is  localized  on  the  outer  atoms,  which  is  confirmed  by  an  experimental 

81  ^  -f* 

study.  The  Ar^  -*  A^  +  Ar  dissociation  energy  is  found  to  be  equal  to 

n  09 

Dg  -  0.20  eV  (Table  4),  close  to  the  experimental  values  of  0.22  eV  ’  and 
8  9 

0.10  eV.  (Throughout  the  manuscript,  we  shall  refer  to  the  theoretical  well 

depths  Dp  as  dissociation  energies,  since  the  zero-point  vibrational  energies 

are  typically  very  small,  i.e.,  -  0.001-0.01  eV.)  Other  calculations  give 

results  a  bit  lower  than  our  dissociation  energy,  namely  0.18  eV,^  0.16^^ 

0.20  eV,^  0.17  eV^  and  0.15  eV.^  The  Zu  ■+  transition  energy  and 

transition  moment  in  the  equilibrium  configuration  are  found  to  be  equal  to 

13  72 

2.26  eV  and  0.8  Dg ,  respectively,  close  to  other  theoretical  values.  ’ 
According  to  our  calculations  for  the  triatomic  Xe^  cluster,  like  for  Ar^,  50X 

of  the  charge  is  concentrated  at  the  central  atom  (Figure  2) ,  in  accordance 

13  73  ’ 

with  other  calculations,  ’  and  the  Xe-Xe  separation  is  equal  to  3.30  A 

(3.47  A  in  Reference  13  and  3.32  A  in  Reference  72).  The  Xe^  -*  Xe^  +  Xe 

dissociation  energy  is  found  to  be  0.20  eV  (Table  5),  a  little  lower  than  the 

87 

experimental  value  of  0.27  eV  but  halfway  between  the  theoretical  values  of 
13  73 

0.12  eV  and  0.36  eV.  For  the  2u  -*  transition  energy  we  obtain  1.87  eV, 
a  bit  higher  than  in  Reference  13  (1.6  eV) . 

The  most  stable  Ar£  cluster,  according  to  our  DIIS  calculation,  is  formed 
bv  a  collinear  Ar^  ion  and  almost  neutral  (q  -  0.002)  Ar  atom  separated  from 

o 

the  nearest  charged  atom  by  3.68  A  (Figure  1).  The  dissociation  energy  of  the 
Art  -*  Ar^  +  Ar  detachment  process  is  found  to  be  equal  to  Dg  -  0.047  eV,  a 
tvpical  energy  of  the  polarization  attraction  (Table  4) .  According  to  the  DIM 
calculation,  ^  the  most  stable  configuration  (dissociation  energy 

D  -  0.043  eV)  has  a  structure  Ar^Ar,  like  in  our  calculation,  but  with  a 

e  i 
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linear  geometry.  In  our  calculation  the  linear  asymmetrical  Ar^Ar  cluster 
presents  a  metastable  structure  with  a  dissociation  energy  of  0.031  eV  only. 
Since  the  differences  in  the  energies  of  the  isomeric  Ar£  structures  are 
small,  most  probably  bevond  the  accuracy  of  the  semiempirical  calculations,  it 
is  really  difficult  to  conclude  which  geometry  is  the  most  stable  one.  The 
more  important  finding  of  the  semiempirical  calculations^ ’ ^ ^  is  the 
asymmetrical  Ar^Ar  structure  of  the  Ar£  cluster.  Due  to  this  structure  the 

Art  excitation  spectrum  is  expected  to  be  similar  to  the  Art  spectrum.  In 

addition  to  the  asymmetrical  Ar^Ar  structures,  we  find  also  a  metastable 
symmetrical  Art  configuration  with  the  geometry  of  an  equilateral  pyramid  with 

o 

a  separation  of  2.836  A  between  every  two  atoms  (F*gure  1,  Table  4).  The 
energy  of  this  fullv  symmetrical  structure  is  only  0.01  eV  above  the  energy  of 
the  most  stable  Ar^  configuration.  In  the  case  of  xenon  clusters  just  the 

symmetrical  equilateral  pyramid  (Xe  detachment  energy  as  large  as  - 

0.203  eV)  is  found  to  be  the  most  stable  configuration  of  the  Xe^  cluster 
(Figure  2).^  The  most  stable  asymmetrical  Xe^Xe  configuration  in  our 
calculation  has  a  bent  structure  and  the  Xe  detachment  energy  of  0.079  eV, 
much  smaller  value  than  that  of  the  symmetrical  pyramid.  This  asymmetrical 
configuration  is  separated  from  the  configuration  of  the  symmetrical  pyramid 
by  a  high  barrier  of  approximately  0.3  eV.  In  contrast  to  our  results,  in 
Reference  73  the  most  stable  Xe£  structure  (Xe  detachment  energy  of  0.159  eV) 
is  a  symmetrical  linear  geometry.  In  our  calculation  the  symmetrical  linear 
Xe^  cluster  is  a  metastabie  configuration  with  a  Xe  detachment  energy  of 
0.069  eV  only.  The  three  Xe*  configurations,  namely  the  symmetrical  pyramid, 
the  symmetrical  linear  structure  and  the  asymmetrical  (bent)  Xe^Xe  structure, 
are  expected  to  have  different  spectrum  features,  so  that  they  can  be 
distinguished  by  experimental  spectral  studies. 

In  the  R*  (n  >  3)  clusters,  the  charge  is  supposed  to  be  localized  at  an 
ionic  core  consisting  of  a  few  R  atoms.  This  ionic  core  is  surrounded  by 
neutral  or  almost  neutral  atoms  which  are  attracted  to  the  charged  core  by 
polarization  and,  to  a  lesser  degree,  by  dispersion  forces  (Tables  4  and  5). 

The  main  problem  which  has  been  intensively  discussed  in  many  experimental  and 
theoretical  works  is  the  size  and  structure  of  the  ionic  cores.  The  simplest 
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model  of  one  atom  ionic  core  is  unrealistic,  as  the  rare  gas  ions  R  are 

chemically  active  and  are  able  to  form  diatomic  molecules  Rj  or  Rn+  clusters 

where  charge  is  shared  by  all  cluster  atoms.  The  strongly-bound  (Dg  -  1  eV) 

diatomic  molecules  provide  the  simplest  realistic  model  for  an  ionic 

core.  ’  Since  these  molecules  can  combine,  as  noted  above,  with  more 

R  atoms  forming  relatively  strongly-bound  (D0  -  0.2  eV)  ionic  trimers  or 

weaklv-bound  R^  clusters,  the  three-  and  four-atomic  ionic  cores  should  be 

considered  as  an  alternative  to  the  two -atomic  ionic  cores.  Although  in  the 

gas-phase  these  clusters  are  more  stable  than  the  R^  molecules,  the  stability 

inside  the  cluster  is  not  obvious,  as  the  energy  of  the  polarization  and 

dispersion  interaction  between  an  ionic  core  and  neutral  atoms  and  between 

neutral  atoms  themselves  is  on  the  same  order  as  the  cluster  dissociation 

energy.  The  delocalization  of  the  charge  on  more  atoms  decreases  the  system 

energy  but  decreases  also  the  attractive  polarization  energy,  which  imposes 

limits  on  the  charge  spread.  This  tendency  of  the  charge  localization  on  a 

small  number  of  atoms  is  reinforced  by  the  relaxation  shift  of  atoms  along  the 

lines  of  attractive  forces  and  the  formation  of  a  stable  ionic  core.  The 

ab  initio  calculation  of  small  (n  -  4-7)  He^  clusters**® ’ ^  shows  that  the  most 

stable  ionic  core  has  the  structure  of  a  linear  symmetrical  trimer,  similar  to 

the  free  He^  cluster.  According  to  semiempirical  calculations  of  Ar*  ^ 

and  Xe  ^  ionic  clusters,  the  ionic  core  may  consist  of  either  three  or 

four  atoms.  As  shown  in  Reference  73,  the  linear  symmetrical  trimer  and 

tetramer  cores  yield  very  similar  energies  for  numerous,  almost  degenerate, 

isomers  of  Xe*  clusters.  Slightly  different  results  are  obtained  for  Ar*  and 

Xe*  clusters  in  Reference  71,  where  the  linear  symmetrical  trimer  core  is 

found  to  be  more  stable  in  small  clusters  (n  <  13)  and  the  linear  symmetrical 

tetramer  core  to  be  more  stable  in  larger  (n  >  13)  ones.  On  the  basis  of  this 

finding,  the  last  model  calculation  of  Ar^  (n  -  3-27)  clusters  was  performed 

93 

for  both  trimer  and  tetramer  (linear  symmetrical)  ionic  cores.  In  our 
calculation,^  in  contrast  to  the  results  of  References  71  and  73,  the  most 
stable  Xe*  clusters  are  formed  by  a  symmetrical  equilateral  Xel  core. 

The  spectroscopic,  photodissociation,  and  ionization  studies  of  R+ 
clusters  are  contradictory  in  their  conclusions  concerning  the  structure  of 
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the  ionic  cores.  In  some  of  this  work  the  conclusion  is  made  that  in  Ar^ 
clusters  the  Ar0  dimer  presents  the  core.  '  ’  In  other  works  the 

experimental  evidence  is  interpreted  in  favor  of  the  trimer  core  in 
Ar+  15,75,81  ^£or  a  crj.tique  of  the  interpretation  of  Reference  81,  see 
Reference  96)  and  in  Xen  clusters.  Recent  studies,  however,  present  more 
evidence  about  the  dependence  of  the  ionic  core  structure  on  the  cluster 
size . 77 ’ ^  According  to  this  work,  the  charge  is  located  on  a  trimer  core 
in  small  clusters  (n  <  13  in  Reference  77)  and  is  spread  out  on  a  tetramer 
core  in  larger  clusters,  in  accordance  with  the  DIM  calculation.^  The  rise 
of  the  charge  delocalization  with  the  increase  of  the  cluster  size  is, 
however,  somewhat  contradictory  with  the  fact  that  in  solid  argon  the  ionic 
charge  is  concentrated  on  a  dimer.  The  conclusions  about  the  structure  of 
the  ionic  core  are  based  mainly  on  the  comparison  of  the  spectroscopic  (and 
some  other)  properties  of  the  ionic  clusters  R^  (n  >  4)  with  the  corresponding 
properties  of  the  R*  molecule  and  the  or  R^  free  ionic  clusters.  It  is 
important  to  note  that  such  comparison  may  lead  to  incorrect  conclusions,  as 
the  spectroscopic  properties  depend  on  the  charge  distribution  not  only  in  the 
ground  state  but  also  in  the  excited  states.  Even  if  the  charge  is 
concentrated  in  the  ground  state  on  a  trimer  core,  for  example,  in  the  excited 
state  the  charge  may  be  partly  transferred  to  other  (neutral  in  the  ground 
state)  atoms,  which  affects  the  excitation  spectrum.  This  point  is  supported 

by  our  semiempirical  calculation  of  the  excitations  of  small  Arn  and  Xe^ 

14 

clusters . 

The  size  and  structure  of  ionic  cores  affect  the  stability  of  the  R* 
clusters.  The  mass-spectroscopic  studies  of  the  ionic  R*  clusters  show  that 
some  of  these  clusters  are  much  more  abundant  than  others,  giving  rise  to  so- 
called  magic  numbers.  These  numbers  found  in  different  works  are  not  always 
the  same.  For  example,  in  He^  clusters  it  is  not  clear  whether  n  -  4  is  the 
magic  number  or  n  -  7  , 10 , 14 , 30 . 99 *  In  argon  and  xenon  clusters  n  -  13  is 

found  as  the  magic  number,  like  in  the  neutral  clusters  with  the  ioosahedral 
svmmetrv .  ^ ' 90 '  This  n  «■  13  argon  cluster  is  not  found,  however,  as  the 

magic  one  in  Reference  103.  In  considering  the  abundance  of  the  ionic 
clusters,  we  have  to  keep  in  mind  not  only  the  stability  of  such  clusters  but 
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also  che  reltive  stability  of  neutral  clusters  where  ionization  produces  the 
ionic  clusters.  After  the  cluster  is  ionized  and  the  initial  ionic  core  is 
formed,  the  cluster  relaxation  to  the  stable  ionic  core  gives  rise  to  the 
energy  excess  which  leads  to  the  loss  of  neutral  atoms  and  dimers. 104  The 
dynamics  of  the  ionized  clusters  fragmentation  affects  significantly  the 
relative  abundance  of  ionic  clusters. 10^ 

VI .  Rare  Gas-Hvdrogen  Ionic  Clusters 

Rare  gas  atoms  can  form  not  only  the  homonuclear  clusters  R*  considered 
in  the  previous  section  but  also  ionic  heteronuclear  clusters  consisting 
either  of  different  rare  gas  atoms  or  rare  gas  atoms  and  molecules  or  valence- 
active  atoms.  The  charge  distribution  and  the  structure  of  heteronuclear 
ionic  clusters  depend  strongly  on  the  relation  between  the  ionization 
potentials  of  che  atoms  involved.  If  the  cluster  is  formed  by  atoms  and 
molecules  with  very  different  ionization  potentials,  the  charge  is  obviously 
localized  on  the  particles  with  lower  ionization  potential.  When  valence- 
active  atoms  are  involved,  the  situation  becomes  more  complicated  because  of 
the  formation  of  valence-bonded  molecules,  such  as  H2 ,  H^,  and  (RH)+  in  the 
case  of  rare  gas -hydrogen  clusters. 

The  rare  gas-hydrogen  systems  (RH2)+  have  been  studied  in  the  context  of 

the  R^  +  H2  and  R  +  chemical  reactions.1’1'0  The  study  of  the  stable 

rare  gas -hydrogen  clusters  has  been  concentrated  mainly  on  the  case  of 

ArH^,114  whose  structure,  by  analogy  with  HeH^,115  is  suggested  to  be  a 

triangle  with  Ar  at  the  vertex.  This  structure  is  supported  by  an  ab  initio 

calculation. 116  The  charge  is  localized  on  the  ion  also  in  (H2^nH3  ionic 

clusters .  1117  ’  11J1  since  the  ionization  potential  of  Hj  (15.34  eV)  is  much 

higher  than  che  H  ionization  potential.  The  H  ionization  potential  is  also 

much  lower  than  the  Ne  ionization  potential,  so  that  in  the  ionic  clusters 

Ne  H4"  the  atom  H  is  expected  to  bear  the  charge,  what  is  confirmed  by  an 
n 

ab  initio  calculation.119  The  rare  gas  atoms  are  also  usually  neutral  (or 
almost  neutral)  in  clusters  which  they  form  with  ionic  molecules,  such  as 
120  or  aromatic  cations . 51 '  53  ’  121  ’ 122  In  all  these  systems  the  charge  is 
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strongly  localized  and  the  particles  are  kept  together  mainly  by  polarization 

forces.  The  heteronuclear  ionic  clusters  with  rare  gas  atoms  has  recently 

123  - 125 

attracted  much  interest.  The  new  (cluster)  class  of  reactions  suggest 

possible  practical  application  of  the  heteronuclear  ionic  clusters. 

Considering  the  structure  of  rare  gas -hydrogen  ionic  clusters  by  the 

14 

semiempirical  DIIS  method,  we  decided  to  look  at  Ar  and  Xe,  since  they  are 
expected  to  form  different  kinds  of  clusters  with  H,  as  the  H  ionization 
potential  (13.6  eV)  lies  between  the  ionization  potentials  of  Ar  and  Xe 
(Table  1) . 

According  to  the  results  of  our  semiempirical  DIIS  calculation,  the 

argon-hydrogen  (Ar^H)^  clusters  are  formed  by  an  (ArH)+  molecule  and  n-1 

neutral  Ar  atoms  which  are  attracted  to  the  ionic  molecule  by  polarization 

forces  (Table  5).  The  (ArH)+  ion  is  a  strongly-bound  molecule  with 
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dissociation  energy  of  4.055  eV  and  equilibrium  distance  of  1.266  A.  ’  ’ 

The  neutral  Ar  atoms  are  located  on  the  H  side  of  the  (ArH)+  molecule,  as  the 
H  atom  bears  larger  (+0.569)  charge  than  the  Ar  atom.  The  Ar(ArH)+  cluster 

o 

has  a  bent  geometry  with  the  neutral  Ar  atom  at  a  distance  of  2.81  A  from  H. 

In  the  Arn(ArH)+  (n  -  2-4)  clusters,  the  neutral  Ar  atoms  are  located 

symmetrically  around  the  (ArH)+  axis.  The  (ArH)+  molecule  is  only  slightly 

affected  by  the  surrounding  neutral  atoms  in  the  cluster.  In  the  sequence  of 

Ar^(ArH)+  clusters,  the  charge  is  shifted  slightly  from  the  Ar  atom  to  the 

H  atom,  but  the  interatomic  distance  is  practically  not  changed.  The  energy 

of  the  Ar  detachment  in  the  Ar^(ArH)+  clusters  is  decreasing  with  k  from 

0.19  eV  for  k  -  3  to  0.09  eV  for  k  -  6  (Table  6).  The  electronic  spectrum  of 

these  was  found  to  be  practically  identical  to  the  free  (ArH)+  spectrum. 

The  results  of  our  semiempirical  calculation  contradict  those  of  an 
127  + 

ab  x.nicio  calculation,  indicating  the  (ArH)  molecule  to  be  strongly 
affected  by  other  Ar  atoms.  The  simplest  argon-hydrogen  cluster  has  a  linear 

o 

ArHAr  asymmetric  structure  with  the  Ar-H  distances  of  1.42  A  (compared  to 

O  o  o 

1.266  A  in  a  free  (ArH)  molecule)  and  1.0  A  (compared  to  2.8  A  in  our 
calculation).  Also,  in  larger  (Ar^H)+  clusters  one  of  the  Ar  atoms  is  bound 
more  strongly  with  the  H  atom  than  others,  but  not  as  strongly  as  in  a  free 
(ArH)+  molecule.  The  ab  initio  results  look  somewhat  strange  for  two  reasons. 
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First,  the  hydrogen  is  a  one-valence  atom  and  cannot  provide  a  valence  bond 
with  two  atoms.  Second,  the  neutral  Ar  atom  polarization  energy  is  on  the 
order  of  0.2  eV,  a  much  smaller  value  than  the  (ArH)+  dissociation  energy. 

The  argon-hydrogen  molecule  ionic  cluster  (Ar^)^"  has  three  asymptotic 
states,  Ar+  +  H0 ,  Ar  +  and  H  +  (ArH)+,  with  dissociation  energies  of 
4.76  eV,  4.95  eV  and  6.215  eV,  respectively.  The  most  stable  cluster  can  be 
formed  by  the  (ArH)+  molecule  and  H  atom.  However,  according  to  our 
calculation,  the  polarization  attraction  of  the  H  atom  to  the  ionic  molecule 
is  so  weak  that  the  H  atom  and  (ArH)+  are  not  attached.  The  lowest 
quasistable  (ArH0)+  cluster  has  a  triangle  geometry  and  is  formed  by  the 
ion  and  almost  neutral  (q  -  +0.02)  Ar  atom  located  symmetrically  with  respect 
to  both  H  atoms.  The  Ar  +  dissociation  energy  is  found  to  be  0.19  eV.  The 
quasistabie  ArH^  cluster  is  separated  from  the  stable  (ArH)+  +  H  state  by  a 
barrier  of,  roughly,  0.75  eV.  The  structures  of  the  quasistabie  Ar^Hj  (k  <  6) 
clusters  are  similar  to  that  of  the  ArH^  cluster,  i.e.,  the  almost  neutral 
Ar  atoms  are  located  in  the  plane  perpendicular  to  the  axis. 

The  xenon- hydrogen  ionic  clusters  (Xe^H)4”,  like  the  argon-hydrogen 

clusters,  are  formed  by  neutral  Xe  atoms  and  the  (XeH)+  ionic  molecule.  The 

(XeH)+  molecule  is  a  strongly-bound  system  with  dissociation  energy  around 

4  eV  and  equilibrium  distances  of  1.6  A.  ^  ^  In  the  Xej{.(XeH)+  clusters,  in 

contrast  to  the  Ar^CArH)4^  clusters,  the  neutral  Xe  atoms  are  located  on  the  Xe 

side  of  the  (XeH)+  molecule,  since  in  this  molecule  the  rare  gas  atom  bears 

most  of  the  charge.  The  Xe(XeH)+  cluster  is  found  to  have  a  linear  XeXeH 

geometry.  The  Xe  detachment  energy  in  the  Xe^(XeH)+  clusters  is  low, 

0.04-0.08  eV  only.  The  (Xe2H)+  clusters  have  been  studied  experimentally  in 
128 

Xe  solids.  The  structures  of  these  clusters  are  found  to  be  different  from 

those  obtained  by  us  for  the  gas  phase. 

The  xenon-hydrogen  molecule  (XeH2)+  ionic  clusters  are  expected  to  be 
formed  by  Xe  +  ion  and  a  neutral  Hj  molecule,  whose  ionization  potential  is 
about  3.3  eV  higher  than  that  of  Xe .  The  DIIS  calculation,  however,  shows  the 
existence  of  a  strongly-bound  linear  (XeHH)+  cluster  with  some  transfer  of  the 
charge  from  the  Xe+  ion  to  the  next  H  atom,  which  bears  a  small  charge  of 

o  o 

■t-0.044  The  Xe-H  distance  in  this  cluster  is  1.79  A,  0.18  A  larger  than  in 
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che  free  (XeH)+  molecule,  and  the  H-H  distance  is  0.81  A,  0.07  A  larger  than 
in  the  free  H2  molecule.  The  energy  of  the  Xe+  +  H2  dissociation  of  the 
cluster  is  large,  0.456  eV.  The  electronic  structure  of  the  (XeHH)+  cluster 
can  be  described  as  a  superposition  of  the  Xe+HH  and  (XeH^H  configurations. 
The  existence  of  the  valence  bound  (XeHH)+  cluster  can  be  detected 
experimentally  by  studying  its  spectrum,  which  is  found  to  be  very  different 
from  that  of  the  H0  and  (XeH)+  molecules.  The  structure  of  the  (XekH2)+ 
clusters  with  more  than  one  Xe  atom  is  quite  different  from  the  (XeHH)+ 
structure.  The  (Xe^Ho)^  clusters  are  formed  by  ionic  h£  clusters  and  a 
neutral  H2  molecule  which  is  bound  to  the  charged  atoms  by  polarization 
forces.  The  Xe£  +  H2  dissociation  energy  of  the  XekH2  clusters  is  about 
0.1  eV.  Due  to  the  interaction  with  the  H2  molecule,  the  Xe£  geometry  in  the 

Xe^Hn  clusters  is  different  from  that  of  the  free  Xe^  clusters. 

k  2  K 

In  the  semiempirical  DIIS  calculation  of  the  rare  gas-halogen  ionic 
clusters,  empirical  or  ab  inicio  potentials  of  the  following  diatomic 
fragments  were  used:  R2 ,  H2 ,  RH,  R2,  H2  and  (RH)  ,  R  ~  Ar,Xe.  These  diatomic 
potentials  were  fitted  by  simple  analytical  expressions  presented  in 
Reference  14. 


VII .  Rare  Gas-Halogen  Systems  with  Electron  Transfer 
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The  rare  gas -halogen  quasistable  molecules  with  electron  transfer  R  X 
and  R^X*  have  attracted  much  interest,  particularly  because  of  their  lasing 
properties. ^9-136  excited  molecules  are  obtained  by  electron  impact  or 

bv  optical  irradiation  of  rare  gas -halogen  mixtures.  The  quasistable 

lecules  R2X'  have  also  been  detected  in  the  condensed  phase. 137-143  jn  a 
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pure  rare  gas  matrix,  for  example  in  Xe  solid  doped  with  Cl2  and  HC1,  near-UV 
irradiation  produces  the  quasistable  triatomic  molecules,  whereas  in  a  mixed 
rare  gas  host,  such  as  Kr  and  Ar  solid  with  a  Xe  admixture,  the  diatomic 
Xe+Cl "  molecules  are  formed. The  R+X*  and  R~X"  molecules  in  solids,  in 


contrast  to  the  gas  phase,  interact  with  other  host  rare  gas  atoms  to  form 
some  systems  which  can  be  described  as  R*X'  (n  >  2)  clusters . ^  More 
than  two  rare  gas  atoms  are  also  involved  in  solids  in  the  excitation  to  the 
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states  with  electron  transfer  from  R  to  X  atom,  which  is  well  demonstrated  by 
excitation  spectra. ^ 

These  rare  gas -halogen  molecules  with  electron  transfer  have  also  been 
studied  theoretical! v .  Ab  inicia  calculations  of  diatomic  molecules  have  been 
performed  for  XeX,  X  **  F, Cl , Br , I ,  in  both  ground  (neutral  configuration)  and 
excited  (ionic  configuration  R+X')  states.  Ab  initio  calculations  of  the 
triatomic  R£X  molecules  have  been  performed  for  R2F,  R  -  Ar.Kr,^^  and 
XejCl.^^  Semiempirical  DIM  calculations  of  K^Cl  and  R2F,  R  -  Ne,Ar,Kr,^^ 
and  XejCl^"^  have  been  performed  for  the  excited  R^X"  (ionic)  configurations 
but  without  taking  into  account  the  coupling  with  the  ground  state  (neutral) 
configuration. 

The  R^X  molecules  and  clusters  in  both  ground  state  (neutral)  and  excited 
(ionic)  configurations  can  be  described  as  systems  with  one  hole  in  their 
closed-shell  atoms.  Such  systems,  as  discussed  in  Section  IV,  can  be 
considered  by  the  semiempirical  DIIS  method  Since  the  RnX  systems  in  the 
excited  states  with  electron  transfer  (ionic  configuration  R^X~)  have  two 
charges,  the  polarization  energy  is  not  ->f  pairwise  character,  which  makes  the 
calculation  more  complicated  (see  Section  III).  It  is  important  to  note  that 
due  to  the  coupling  between  the  neutral  and  ionic  configurations,  some  small 
electron  transfer  from  R  to  X  takes  place  also  in  the  neutral  state  (see 
Section  II),  with  the  electron  transfer  in  the  excited  R^X"  states  not  being 
lOOi.  The  DIIS  method  takes  into  account  just  this  coupling,  opposite  to  the 


DIM  method  in  the  version  used  for  the  calculation  of  the  R2X 
molecules .  ^  The  DIIS  method  has  been  applied  by  us  for  the 

o  rn  IQ 

semiempirical  calculations  of  Xe^Cl  ’  and  Xe^HCl  clusters  in  the  gas 

,  8,39  .  39,60 

phase  and  solids. 

The  DIIS  calculation  of  the  triatomic  quasistable  mole  ule  Xe^Cl 


provides  potential  energy  surfaces  similar  to  those  of  the  DIM 

calculation  Unlike  DIM,  the  DIIS  calculation  provides  also  the  dipole 


transition  moments  between  excited  and  ground  states.  The  lowest  quasistable 

2 

excited  state,  4  F.  according  to  the  DIIS  calculation.  has  a  configuration  of 


an  isosceles  triangle  with  the  Xe-Cl-Xe  angle  of  60.4°  and  Xe-Cl  distance  of 

°  °  + 

3.23  A.  0 . 16  A  larger  than  in  the  quasistable  diatomic  Xe  Cl  excimer.  The 
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“  + 

Xe-Xe  distance  is  3.25  A,  almost  the  same  as  in  the  ionic  Xe2  molecule. 

Qualitatively,  the  excited  trimer  can  be  described  as  a  Xe^  molecule  and  C1‘ 

ion  which  are  attracted  to  each  other  by  a  Coulombic  force.  The  coupling 

between  the  ionic  and  neutral  configuration  decreases  slightly  the  Xe£  and  Cl" 

2 

charge,  which  in  the  4  T  state  is  found  to  be  0.964  and  -0.964,  respectively. 

2 

The  coupling  decr-ases  the  energy  of  the  4  T  quasistable  trimer  by  0.07  eV  and 

increases  the  ground- state  energy  in  the  geometry  of  the  quasistable  trimer  by 

2 

0.1  eV.  The  energy  of  the  rransition  from  the  quasistable  4  T  state  to  the 

2 

1  T  ground  state  is  found  to  equal  2.56  eV  (435  nm) ,  which  is  close  to  the 

2  o  136  22 

center  of  the  4^F  -»  l'T  emission  band.  The  A  T  -*■  IT  transition  moment  has 

a  sharp  minimum  in  the  Xe^Cl"  equilibrium  geometry. 

A  quasistable  molecule  with  rare  gas-halogen  electron  transfer  is  found 

2  2 

also  in  the  n-svmmetry  4  F^  state.  Like  in  the  E-symmetry  4  T  state,  the  II- 

svmmetrv  quasistable  molecule  has  the  geometry  of  an  isosoeles  triangle,  but 

with  a  much  larger  Xe-Cl-Xe  angle  of  85°  .  The  Xe-Cl  distance  in  this  trimer 

is  3.13  A,  shorter  than  in  the  4  F  state,  whereas  the  Xe-Xe  distance  is  much 

2  2 

larger  than  in  the  4  T  state,  so  that  the  4  Tn  quasistable  molecule  cannot  be 

considered  as  formed  by  a  Cl"  ion  and  Xe^  ionic  molecule.  The  energy  of  the 

2 

quasistable  4  molecule  is  0.8  eV  higher  than  that  of  the  lowest  quasistable 

2 

4  T  molecule. 

There  are  *"vo  possible  models  of  the  structure  of  quasistable  Xe^Cl" 

(n  >  2)  clusters.  First,  they  can  be  formed  by  the  Xe^Cl  trimer  and  neutral 

Xe  atoms,  which  are  attracted  to  the  trimer  by  polarization  forces.  The 

polarization  forces  in  this  case  are  expected  to  be  weaker  than  in  • he  case  of 

ionic  clusters  (see  Sections  V  and  VI) ,  since  the  Xe^Cl"  trimer  is  not  a 

charged  particle  but  a  particle  with  a  dipole  moment.  The  second  possibility 

is  the  formation  of  the  Xe+Cl  cluster  by  the  Cl  ion  and  the  ionic  Xe+ 

n  n 

cluster,  especially  the  relatively  stable  Xe^  cluster  which,  in  a  free  state, 
has  most  probably  a  symmetrical  linear  structure  (see  Section  V).  According 
to  the  results  of  our  calculation  for  n  -  3,  the  first  possibility  is 
realized,  and  the  most  stable  Xe^Cl  cluster  has  the  XeXe^Cl"  structure  with 
the  neutral  Xe  atom  located  symmetrically  above  the  triangle  Xe^Cl  plane  at  a 

a  o 

distance  of  3.52  A  from  the  Cl  atom  and  4.41  A  from  two  charged  Xe  atoms. 
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This  Xe  atom  is  noc  completely  neutral,  but  the  charge  it  bears  is  negligibly 
small,  +0.005  only.  The  energy  of  the  neutral  Xe  atom  detachment  is 
relatively  large,  0.11  eV.  The  XeXejCl  vertical  transition  energy  to  the 
ground  state  is  2.48  eV,  close  to  the  XejCl'transition  energy  of  2.56  eV, 
which  makes  it  difficult  to  detect  the  XeXe^Cl"  clusters  by  their  emission 
energy.  The  less  stable  cluster  (Xe  detachment  energy  of  0.06  eV)  can  be 
described  as  formed  by  the  Xe^  cluster  and  Cl"  ion;  however,  the  Xe-Xe 
distances  in  this  Xe^Cl  cluster  (3.58  A)  are  much  larger  than  in  the  free  Xe^ 
cluster  (3.38  A),  and  consequently  the  three  Xe  atoms  do  not  form  a  valence 

bond  system.  The  Cl  ion  is  located  symmetrically  with  a  distance  from  the 

°  , 
central  Xe  atom  of  3.33  A.  Like  in  the  n  —  3  case,  the  most  stable  Xe  Cl" 

n 

(n  >  3)  clusters  are  formed  by  the  Xe^Cl"  trimer  and  almost  neutral  n-2  Xe 
atoms.  The  structure  of  the  n  -  12  cluster,  for  example,  is  Xe^Xe^Cl" .  The 
Xe^Cl"  trimer  in  this  cluster  is  a  little  bit  deformed,  and  its  Xe-Xe  and  Xe- 

o 

Cl  distances  (equal  one  to  another)  are  larger  by  0.04-0.06  A  than  in  the  free 
XejCl  cluster.  This  deformation  of  the  Xe^Cl  trimer  is  mainly  due  to  some 
loss  of  the  positive  charge,  which  is  transferred  to  the  ten  so-called  neutral 
atoms  which  bear  together  4.3X  of  the  positive  charge. 

The  Xe^Cl  system  is  of  great  interest  in  regard  to  Xe  solids  doped  by 
Cl-containing  molecules  .  ^  These  molecules,  usually  C^  or  HC1,  are 

dissociated  by  irradiation  producing  free  Cl  atoms  which  replace  Xe  atoms  and 
form  a  CLXe-^  cage  with  the  Cl  atom  in  the  middle.  The  photon  excitation  of 
this  cage  leads  to  electron  transfer  from  the  Xe  atoms  to  the  Cl  atom  and  the 
formation  of  the  spherically  symmetrical  or  almost  symmetrical  Xe^Cl"  complex 
which,  in  accordance  with  the  results  obtained  for  the  Xe^Cl  cluster,  is 
unstable  and  tends  to  form  a  Xe^gXe^Cl"  complex.  According  to  the  DIIS 
calculation  of  this  complex  in  its  equilibrium  geometry,  the  Cl"  atom  is 

o 

s  .ifted  from  the  center  of  the  Xe  matrix  cage  by  1.12  A  forming  with  two  Xe 
atoms  the  Xe^Cl  trimer.  The  trimer  has  the  shape  of  an  isosoeles  triangle 

o 

with  the  quasistable  Xe-Cl-Xe  angle  as  60.2°,  each  Xe-Cl  distance  as  3.30  A, 

o 

and  the  Xe-Xe  distance  as  3.31  A.  These  distances  are  larger  than  in  the  free 
+  •  ° 

Xe9Cl  cluster  (3.22  A  and  3.25  A,  respectively)  but  are  close  to  the 

+  -  ° 

interatomic  distances  in  the  Xe^gXe^Cl  cluster  (3.29  A).  The  ionic  charge 
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of  the  Cl'  ion  is  -0.97,  almost  the  same  as  in  the  gas-phase  XejCl"  molecule 
(-0.964).  The  positive  charge  is  partly  delocalized,  so  that  the  common 
charge  of  the  two  Xe  atoms  of  the  trimer  is  +0.932,  whereas  +0.038  of  the 
positive  charge  is  distributed  among  the  other  ten  atoms  of  the  cage.  One  of 
these  atoms  which  carries  most  of  this  charge,  +0.021,  is  shifted  further  from 

o 

the  Cl'  ion  by  0.1  A.  Nine  other  Xe  cage  atoms  are  located  at  the  same  points 
as  in  the  ground  (neutral)  state.  The  transition  from  the  quasistable  Xe£  C1‘ 
configuration  to  the  ground  state  has  a  large  transition  moment  of  2.41  Dg. 

The  vertical  transition  energy  is  2.95  eV,  0.21  eV  smaller  than  the  gas-phase 
Xe^2Cl"  transition  energy.  The  experimental  value  for  the  shift  of  the 
Xe^Cl'  -  Xe2Cl  transition  energy  in  the  solid,  compared  to  the  gas  phase,  is 

C  Q 

larger,  namely  0.38  eV.  The  calculation  of  the  potential  curves  for  the 
inner  motion  of  the  atoms  inside  the  Xe^Cl ’  trimer  allows  us  to  estimate 
vibrational  excitations,  which  are  found  to  be  small,  0.034  eV  and  0.012  eV. 
The  vibrational  motion  broadens  by  some  0.28  eV  the  energy  of  the  transition 
to  the  repulsive  ground  state  and  affects  significantly  the  shape  of  the 
emission  spectrum. 

The  rare  gas  atoms  are  known  to  form  van  der  Waals  (vdW)  clusters  with 

HC1  molecule.37,1^9'150  These  clusters  can  be  excited,  in  particular  in  the 

rare  gas  solids,  to  states  with  electron  transfer  forming  quasistable 

complexes  similar  to  the  R^X'  complexes  considered  above.  In  order  to  study 

the  possibility  of  the  formation  of  such  quasistable  complexes,  we  have 

39 

performed  DIIS  calculations  for  XenClH  systems.  According  to  our  results, 

strongly  bound  XeClH  quasistable  molecules  with  electron  transfer  are  formed 

in  several  excited  states  and  surprisingly  even  in  the  ground  state 

(Figure  3).  The  ground  state  molecule  is  found  to  have  a  collinear  (HXe)+Cl 

configuration  of  Z  symmetry  with  a  very  large  dipole  moment  of  15  Dg .  The 

energy  of  this  molecule  is  low,  -0,39  eV,  only  4.24  eV  higher  than  the  energy 

of  the  vdW  complex.  The  (HXe)+Cl'  molecule  is  separated  from  the  XeHCl  vdW 

complex  by  a  barrier  of  0.41  eV  and  consequently  can  be  detected  at  low 

temperatures.  In  the  first  excited  state  the  linear  Z  symmetry  quasistable 

molecule  is  also  formed,  but  its  structure  is  different,  namely  (XeH)+Cl  . 

^  » 

This  molecule  is  separated  from  the  excited  Xe(HCl)  vdW  complex  by  a  barrier 
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of  0,18  eV.  The  (XeH)+Cl’  molecule  emits  in  the  near-UV  (232  pm).  In  the 
higher  excited  states  the  quasistable  molecules  have  a  structure  of  a  Xe+Cl' 
molecule  and  almost  neutral  H  atom  bound  to  the  Xe+Cl’  molecule  by 
polarization  forces.  The  energy  of  the  H  atom  detachment  from  this  molecule 
is  relatively  large,  0.12-0.18  eV.  The  Xe+C1*H  transition  energy  to  the 
ground  state  lies  within  the  limits  of  7. 3-7.5  eV.  In  one  of  the  excited 
states  a  (XeH)+Cl  molecule  is  formed  whose  transition  energy  to  the  ground 
state  is  very  low,  0.8  eV.  It  follows  that  the  quasistable  XeHCl  molecule 
provides  emission  both  in  the  near  infrared  and  in  the  UV  regions  and  presents 
an  interesting  object  for  spectroscopic  study.  The  quasistable  XeHCl 

molecules  with  electron  transfer  can  attract  more  Xe  atoms  forming  quasistable 

39 

clusters . 

In  the  Xe  solids  the  HCl  molecule  is  photodissociated  as  a  result  of 
valence  excitation  to  the  n  state.  It  is  not  immediately  obvious  that  after 
dissociation  the  H  atom  is  moving  to  other  cages,  thus  preventing  the  H-Cl 
association.  However,  according  to  our  calculation  the  H  atom,  after  its 
detachment  from  Cl,  can  surmount  a  barrier  which  separates  adjacent  cages  and 
escape  to  the  bulk,  leaving  the  free  Cl  atom  inside  a  cage,  in  accordance  with 
experiments  .  ^ 

The  excited  rare  gas -halogen  systems  with  electron  transfer  are  of 
interest  not  only  in  connection  with  the  quasistable  systems  considered  above, 
but  also  from  the  point  of  view  of  the  electron  excitation  of  the  ground-state 
rare  gas -halogen  system,  such  as  vdW  molecules  and  clusters  or  the  ClXe^ 
cages  in  the  solid  matrix.  The  problem  of  the  rare  gas-halogen  systems 
excitation  is  of  great  interest,  particularly  in  solids  where  the  excited 
species  interact  with  the  matrix  atoms  .  ^  jn  Xe  matrix  the  Cl  atom 

has  some  freedom  of  motion  so  that  the  excitation  may  take  place  in  different 

49 

localizations  of  the  Cl  atom.  Since  the  energy  and  the  transition  moment  of 
the  excitation  is  affected  by  the  Cl-Xe  atoms  interaction,  which  depends  on 
the  Cl  atom  location,  the  quantum  and  thermal  motion  of  Cl  determines  to  some 
extent  the  excitation  spectrum.  The  DIIS  calculation  performed  for  the 
ground-state  ClXe^  system  in  the  Xe  matrix  shows  that  the  potential  energy 
surface  of  the  Cl  atom  inside  the  cage  has  a  broad  region  of  almost  constant 
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potential,  with  variations  within  the  range  of  1  meV.  The  potential  energy 
surface  of  the  Cl  motion  has  twelve  very  shallow  potential  wells  which  are 

o 

shifted  by  0.5 -0.8  A  from  the  cage  center  toward  the  Xe  atoms  (the  Xe  distance 

o 

from  the  center  is  4.33  A).  The  barriers  between  these  wells  are  so  low  that 
even  quantum  effects  allow  the  atom  to  move  inside  a  relatively  wide  region  of 
the  cage.  Thermal  effects  obviously  make  the  motion  stronger  and  increase  the 
cage  region  accessible  for  the  atom.  Taking  into  account  both  quantum  and 
thermal  motion,  it  is  possible  to  determine  the  probability  of  the  Cl  atom 
location  as  a  function  of  its  location  inside  the  cage.  Calculating  also  the 
excited  states  energy  and  the  transition  moments  for  any  Cl  atom  location,  it 
is  possible  to  determine  the  excitation  spectrum.  According  to  the  results  of 
this  calculation,  the  10  K  and  80  K  spectra  do  not  differ  much  from  other, 
Leading  us  to  the  conclusion  that  the  thermal  motion  does  not  affect  photon 
absorption  to  any  significant  extent.  The  main  contribution  to  the  absorption 
provides  the  transition  to  the  second  excited  state.  The  spectra  demonstrate 
two  wide  maxima  at  320-340  nm  and  270-290  nm  which  are  separated  by  a  deep 
minimum  (Figure  4) .  The  first  of  these  maxima  is  found  to  be  higher  than  the 
second  one.  The  maxima  and  the  lower  limit  of  the  absorption  are  shifted 
slightly  to  smaller  energies  when  the  temperature  increases.  Our  calculated 
excitation  spectrum  is  similar  to  the  experimental  spectrum,  which  also 
demonstrates  two  broad  maxima,  roughly  in  the  same  wavelength  intervals.  In 
accordance  with  the  theoretical  prediction,  the  low-energy  maximum  in  the 
experimental  spectrum  is  larger  than  the  high-energy  maximum. ^ 

VIII.  Rvdberg-Excited  Rare  Gas  Clusters 

In  the  Rydberg-excited  rare  gas  atoms,  the  excited  electron  distance  from 
the  atom  center  is  larger  on  average  than  the  interatomic  distances  between 
neighboring  atoms,  which  suggests  that  the  interatomic  R  -R  potential  is 
determined  mainly  by  the  interaction  between  the  ionic  core  R+  and  neutral 
R  atom.  This  suggestion  is  supported  by  the  existence  of  metastable  Rydberg 
excited  rare  gas  dimers  R^ ,  resembling  ionic  dimers  Rj.  "  The  Rj 

potential  energy  curves,  found  both  experimentally^'7  *  and  by  model^^ 
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and  ab  initio1^ calculations  for  low-lying  states,  demonstrate, 

however,  significantly  weaker  binding  than  in  the  ionic  Rj.  The  relatively 

weak  binding  between  the  excited  and  neutral  atoms  is  attributed  mainly  to  the 

repulsive  forces  arising  between  the  excited  electron  and  the  electrons  of 

another  atom. 172-173  s^nce  the  Rydberg-excited  R*  atom  chemical  activity  is 

weaker  than  that  of  the  R+  ion,  the  ability  of  the  R*  atom  to  form  Rn  (n  >  2) 

clusters,  similar  to  the  R*  clusters,  poses  a  question.  According  to  our 

knowledge,  there  is  not  any  experimental  evidence  about  the  existence  of 

metastable  Rydberg  clusters.  Experimental  studies  of  Rn  Rydberg  states  have 

been  restricted  to  the  problem  of  the  excitation  of  vdW  R2  dimers1^  1'?6  and 

Rn  clusters1'*6' 179  to  autoionizing  Rydberg  states.  For  the  case  of  Ar  dimers 

and  clusters,  the  autoionizing  states  are  formed  by  the  excitation  of  inner- 

176  179 

shell  electrons,  such  as  3s  or  2p.  The  3s  -»  np  spectrum  of  the  Ar^ 

cluster  is  found  to  be  similar  to  that  of  Arj  dimer  but  with  wider  resonance 
lines  . 176  This  finding  cannot  be  extended  automatically  to  the  case  of  the 
excitation  of  the  outer  3p  electrons,  which  are  affected  more  strongly  by 
neighboring  atoms.  The  direct  study  of  Rydberg  excitations  by  the  measurement 
of  absorption  spectra  were  performed  for  dimers  in  Reference  162.  Much  more 
is  known  about  the  Rydberg  excited  states  in  rare  gas  solids.16  The 
information  concerning  the  excitations  in  solids  can  be  useful  for  the  study 
of  large  R  clusters,  but  not  of  the  small  ones.  Any  ab  initio  or 
semiempirical  calculations  of  the  Rydberg  excited  Rn  (n  >  2)  clusters  are 
unknown  to  us . 

The  Rydberg  excited  Rn  systems,  by  analogy  with  the  Rn  clusters 
(Section  V)  or  R^X  systems  (Section  VII),  are  made  by  n  closed-shell  atoms 
with  one  hole  which,  generally  speaking,  is  moving  between  atoms.  Such 
systems,  as  stated  in  Section  IV,  can  be  treated  by  the  semiempirical  DIIS 
method.  Since  the  time  of  the  excited  electron  motion  around  the  ionic  core 
ls  much  shorter  than  the  time  of  the  hole  jump  from  one  atom  to  another,  the 
excited  electron  is  adiabatically  linked  to  the  hole.  It  follows  that 
applying  the  DIIS  method  to  the  excited  clusters  we  have  to  consider,  instead 
of  the  hole  motion,  the  motion  of  the  hole-electron  exciton.  Neglecting  the 
d ipo 1 e - d ipo le  interaction,  which  is  expected  to  be  weak  in  Rydberg-excited 
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systems,  we  are  left  with  pairwise  interactions  which  make  it  possible  to  use 

the  simplified  DIIS  Eq .  (IV. 3).  The  DIIS  method  has  been  applied  by  us  for 

.  ,  ,  _  *  _  180 
the  calculation  of  Ar^  systems. 

In  order  to  calculate  the  electronic  structure  of  Arn  clusters,  one  needs 
as  input  the  Ar2  and  Ar.,  potentials.  The  Ar2  vdW  potential,  used  in  the 
calculation,  fits  the  empirical  equilibrium  distance  and  dissociation 
energy^^  and  the  ab  inicio  potential^"^  at  small  interatomic  separations. 

The  Rydberg-excited  Ar0  potentials  are  taken  from  the  recent  ah  initio 
calculation. 169  The  Ar£  potential  curves  asymptotic  to  the  Ar*(3p54s^)  state 
are  calculated  in  Reference  169  both  without  and  with  spin-orbit  (SO) 
coupling,  whereas  the  Ar2  potential  curves  asymptotic  to  the  Ar  (3p  4p  )  state 
are  calculated  without  SO  coupling  only.  As  we  restrict  ourself  to  the 
Rvdberg  excited  atomic  states  3p54s1  and  3p54p^,  we  shall  denote  further  these 
states  simply  by  the  excited  electron  orbitals  4s  and  4p. 

According  to  the  results  of  the  DIIS  calculation  presented  in  Table  7, 
the  4s  excited  atoms  in  the  states  *2*  and  32*  without  SO  coupling  and  C(0*) 

■jlf 

with  SO  coupling  really  form  metastable  Ar^  clusters.  All  of  these  clusters 

have  a  symmetrical  linear  geometry,  like  the  Ar^  clusters  (Section  IV),  with 

the  Ar-Ar  distances  slightly  larger  than  in  the  corresponding  states  of  the 

Ar2  dimer.  However,  in  contrast  to  the  Ar^  cluster  with  a  dissociation  energy 

Dg  of  0.2  eV,  the  dissociation  energies  of  these  metastable  Ar^  clusters  are 

small,  ranging  in  the  interval  0.01-0.04  eV.  In  the  state  asymptotic  to 

Ar*  (4s ) ,  in  addition  to  the  collinear  metastable  cluster,  a  less-stable 

cluster  with  the  geometry  of  an  equilateral  triangle  is  formed.  The  Rydberg 

3  3 

4s  excited  atom  in  orther  states  with  SO  coupling  (  and  P2)  do  not  form 
any  metastable  clusters,  in  spite  of  the  fact  that  they  form  relatively 
strongly-bound  metastable  dimers  (Table  7).  The  Ar*  clusters  formed  by 
Ar*(4s)  atoms  attract  more  atoms  but  with  so  small  dissociation  energy 
(0.001-0.002  eV)  that  the  very  existence  of  the  metastable  Arn  (n  >  3) 
clusters  asymptotic  to  Ar(4s)  is  under  question. 

The  behavior  of  the  Ar*  atom  excited  to  the  Rydberg  4p  state  shows  much 
more  similarity  to  the  behavior  of  the  Ar+  ion,  which  is  easily  explained  by 
the  more  diffuse  character  of  the  4p  orbital  compared  to  the  4s.  According  to 
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the  DIIS  calculation,  the  dissociation  energy  Dg  of  the  metastable  Ar^  cluster 

3 

in  the  state  asymptotic  to  Ar*  Dg(4p)  is  0,23  eV  (Table  7)  and  exceeds  even 
the  A^  dissociation  energy  of  0.2  eV  (Section  V).  This  Ar^  cluster  attracts 

if 

more  Ar  atoms,  forming  larger  metastable  Ar^  clusters,  like  in  the  case  of 
ionic  clusters. 

Calculating  the  energy  of  the  Rydberg  excited  Ar*  clusters  in  the 
geometry  of  the  ground  state  (vdW)  Arn  clusters,  we  can  determine  the  Arn 
vertical  excitation  energies  (for  the  vdW  clusters  geometry  see  Section  II). 
This  energy  is  determined  as  the  difference  between  the  energies  of  the 
excited  and  ground  states  in  the  same  geometry.  According  to  the  results  of 
the  calculation,  the  energy  of  the  cluster  dissociation  to  separated  atoms  in 
the  Rydberg  excited  states  is  on  the  same  order  as  in  the  ground  states,  so 
that  both  of  these  dissociation  energies  are  of  importance  for  the 
determination  of  the  excitation  energy.  The  dissociation  energy  of  the  ground 
state  (vdW)  clusters  Arn  increases  monotonically  with  n  (Section  II),  coming 
up  to  0.49  eV  for  the  "magic"  Ar-^  cluster.  The  dependence  of  the  excited 
states  dissociation  energy  on  n  is  more  varied.  In  the  states  asymptotic  to 
the  Ar*(4s)  states  and  (without  SO  coupling)  and  ^P^  (with  SO  coupling), 
the  energy  of  the  excited  states  decreases  monotonically  with  the  number  of 
atoms  n,  almost  at  the  same  rate  as  the  vdW  clusters,  which  makes  the 
excitation  energy  dependence  on  n  very  weak.  For  example,  in  the  case  of  the 
cluster  excited  state  asymptotic  to  the  atomic  ^P^  state,  the  excitation 
energy  decreases  slightly  with  n,  up  to  n  -  4-6,  and  increases  for  n  >  6,  but 
slowly  so  that  the  Ar-^  excitation  energy  is  only  0.14  eV  higher  than  the 
atomic  Ar  excitation  energy.  The  situation  is  different  in  the  triplet  states 

3  3 

with  SO  coupling  asymptotic  to  the  Ar  (4s)  P^  and  Pj  atomic  states,  where 
the  variations  in  the  excitation  state  energy  dependence  on  n  are  small.  As  a 
result,  the  excitation  energy  is  increasing  monotonically  with  the  number  of 
atoms  n,  due  to  the  decrease  in  the  energy  of  the  vdW  clusters.  For  example, 

the  Ar^  excitation  energies  to  the  cluster  states  asymptotic  to  the  atomic 

3  3 

and  P.i  states  are  larger  than  the  corresponding  Ar  atomic  excitation 
energies  by  0.34  eV  and  0.29  eV,  respectively.  Unfortunately,  there  are  not 
anv  experimental  aata  on  the  vdW  Arn  clusters  3p  -»  4s  excitations. 
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The  increase,  if  weak,  of  che  vdW  cluster  excitation  energy  with  the 

cluster  size  is  in  contrast  to  the  behavior  of  the  cluster  ionization 

potential.  Due  to  the  polarization  energy  and  hole  delocalization,  the  energy 

of  the  ionic  cluster  (in  the  geometry  of  the  vdW  cluster)  decreases  with 

cluster  size  much  more  strongly  than  the  vdW  cluster  energy,  causing 
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relatively  strong  decrease  in  the  ionization  potential.  For  example,  we 

have  found  the  ionization  potential  of  Ar^j  to  equal  14.46  eV,  smaller  by  1.41 
eV  than  the  atom  Ar  ionization  potential  (15.51  eV) .  It  is  interesting  to 

«i^* 

note  that  in  the  lowest  excited  Ar^  states,  the  excited  electron  is 
delocalized  between  the  12  outer  Ar  atoms,  leaving  the  central  atom  almost 
unexcited.  In  contrast  to  this,  in  the  ionic  Ar^  cluster  (in  the  geometry  of 
the  vriW  Ar, i  cluster),  the  charge  is  concentrated  on  the  central  atom. 

J*  J 

IX.  Concluding  Remarks 

The  rare  gas  atom  R  becomes  chemically  active  when  it  looses  an  outer 
shell  electron.  In  clusters,  the  electron  removed  from  the  atomic  shell 
is  either  eliminated  from  the  cluster  or  left  inside  it.  In  the  first  case, 
an  ionic  cluster  is  formed,  for  example  R*.  The  second  case  can  be  realized 
in  any  cluster  by  Rydberg  excitation  (Rn_iR+e  complex)  and  in  a  halogen 
X-containing  cluster  excited  by  R  -*  X  electron  transfer  (RnlR+X‘  complex). 

Due  to  the  chemical  activity  of  the  R+  ion  or  ionic  core,  R+-R  valence 
attraction  arises  which  leads  to  the  formation  of  a  strongly -bound  R^  or  R* 

(n  >  2)  fragment.  The  ionic  fragment  attracts  other  atoms  by  polarization 
forces.  These  two  main  interactions,  namely  the  valence  binding  resulting 
from  the  positive  charge  delocalization  and  the  polarization  attraction,  are 
taken  into  account  by  our  semiempirical  diatomics-in-ionic-systems  (DIIS) 
method,  which  has  been  used  to  treat  the  rare  gas  clusters  and  solids 
containing  charged  atoms. 

The  structure  and  properties  of  systems  formed  by  two  rare  gas  atoms, 
such  as  R~,  R*  and  Rtx" ,  are  well  known.  We  have  less  reliable  knowledge, 

*.  t- 
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unfortunately,  about  the  systems  R  ,  Rr  and  RnX  with  three  or  more  rare  gas 
From  experimental  and  theoretical  studies,  it  is  known  that  the  ionic 


atoms . 
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clusters  R+,  exist,  but  the  structure  of  these  clusters  is  still  not 
n 

unequivocally  clear.  However,  most  calculations,  including  our  DIIS 
calculation,  have  shown  the  smallest  ionic  clusters,  R^",  to  be  of  a 
symmetrical  linear  geometry.  In  other  calculations  the  geometry  is  found 
to  be  asymmetrical.  The  experimental  evidences  concerning  the  structure 
are  also  contradictory.  The  explanation  of  some  ambiguity  in  the  R^  structure 
lies,  probably,  in  the  features  of  its  potential  energy  surfaces  which  are 
almost  flat  over  wide  intervals  of  interatomic  distances,  so  that  the  thermal 
motion  easily  shifts  the  cluster  geometry  far  away  from  the  equilibrium 
configuration.  The  same  holds  to  an  even  in  a  greater  extent  for  the  R* 
clusters  with  more  than  three  atoms,  which  have  several  isomers  with  almost 
the  same  energy.  The  situation  is  similar  in  the  case  of  clusters  with  R  -*•  X 
electron  transfer,  such  as  Rn_2R2G^'  clusters.  In  the  Rydberg  excited 
clusters,  the  valence  and  polarization  forces  are  less  strong  than  in  the 

C  ■» 

ionic  clusters,  at  leas-  .  the  lowest  np  (n+l)si  excited  state.  According  to 
our  calculations,  th_  ory  existence  of  such  excited  Rn  (n  >  2)  clusters  is 
under  question.  In  higher  excited  states  the  R  -R  interactions  look  more  like 
the  R+-R  interactions ,  so  that  more  similarity  is  expected  in  the  structure  of 
highly-excited  Rn  and  ionic  Rn  clusters. 

Ir  ionic  clusters  R*  with  more  than  four  atoms,  the  charge  is 
concentrated  on  a  smaller  number  of  atoms,  which  form  an  ioni^  core  so  that 
the  structure  of  a  cluster  is  The  size  and  geometry  of  this  ionic 

core  depends,  most  likely,  on  the  cluster  size  n.  The  results  of  both 
.heoretical  and  experimental  studies  concerning  the  size  and  the  geometry  of 
the  ionic  cores  are  contradictory.  According  to  our  calculation,  in  the  Xe* 
clusters  the  ionic  core  has  the  configuration  of  a  fully  symmetrical 
equilateral  XeT  pyramid.  In  the  rare  gas  solids  doped  by  halogen  atoms  in  the 
excited  states  with  R  -»  X  electron  transfer,  the  R^Cl’  molecule  plays  the  role 


of  the  ionic  core. 
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Table 


o 

Atomic  radii  r  (in  A),  ionization  pot'  _>.s  I  (in  eV)  and 
polarizabilities  a  (in  A  ) . 


Atom 

r 

I 

a 

He 

1.33 

24.58 

0.205 

Ne 

1.60 

21.56 

0.395 

Ar 

1.92 

15.755 

1.640 

Kr 

1.98 

14.00 

2.487 

Xe 

2.18 

12.13 

4.017 
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Experimental  interatomic  equilibrium  distances  R  (in  A)  and 
dissociation  energies  DQ  (in  eV)  of  van  der  Waals  R2  molecules. 
(These  are  taken  from  Reference  29.) 


Molecule : 

He„ 

Ne? 

Ar2 

Kr7 

Xe2 

R 

2.96 

3.11 

3.76 

4.01 

4.37 

D0 

0.0024 

0.0095 

0.031 

0.052 

0.085 
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Table  J.  Results  of  model  calculations  for  the  equilibrium  distances  between 

o 

neighboring  atoms  R  (in  A),  dissociation  energies  Dg  and  atom 
detachment  energies  (in  eV)  of  van  der  Walls  Arn  clusters. 


n : 

2 

3 

4 

6 

9 

13 

14 

R 

3.75 

3.75 

3.74 

3.725 

3.72 

3.71 

3.71 

D 

e 

0.012 

0.036 

0.072 

0.153 

0.265 

0.489 

0.544 

da 

0.012 

0.024 

0.036 

0.051 

0.065 

0.070 

0.055 
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Table  4.  Structure  of  the  ionic  clusters  Ar^  and  ionic  molecule  Ar^.  The  distances 

© 

(R)  are  in  A,  and  the  ground-state  and  dissociation  energies  E  and  Dg 


re 

spectively,  in  eV. 

(These 

are  taken 

from 

Reference 

14.) 

N 

Structure  ^ 

R  2> 

r3> 

4) 

q 

E 

Dissociation 

Products 

De 

I 

Ar£ 

2.48 

-1.240 

Ar+  +  Ar 

1.240 

II 

Ar^ 

linear 

2.59 

-1.443 

Ar£  +  Ar 

0.203 

III 

< 

^Ar3  ^  z  ^Ar ^xz 

2.59 

3  .  68 

0.002 

-1.490 

Ar^  +  Ar 

0.047 

IV 

Art 

regular  pyramid 

2.836 

-1.480 

Ar^  +  2At 

0.240 

\j 

Art 

-4 

linear 

2.59 

3  .  3 

0.01 

-1.474 

Atj  +  Ar 

0.031 

VI 

Ar5 

^ Ar3 ) z ( Ar2 } xy 

2.59 

3.68 

0 

-1.540 

(Ar4}  HI  +  Ar 

0.050 

VII 

Ar5 

(Ar3)Z(Ar2)xz 

2.59 

3.61 

0.002 

-1.540 

(Art)m  +  Ar 

0.050 

VIII 

Ar5 

(Ar^z(Ar2)xz 

2.59 

3.68 

0.003 

-1.537 

+  Ar 

0.047 

IX 

Ar6 

(Ar3)2(Ar2)xy(Ar)x 

2.59 

3.61 

0.001 

-1.601 

(Ar^)^^  +  Ar 

0.061 

^The  subscripts  z,  xz ,  etc.  denote  the  location  along  an  axis  or  in  a  plane 
(see  Figure  1). 

2 ) 

Distance  between  adjacent  charged  atoms. 

3 ) 

Distance  from  a  neutral  or  weakly-charged  atom  to  the  nearest  charged  atom. 
^Common  charge  of  all  weakly-charged  atoms. 
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Table  5.  Structures  of  the  ionic  clusters  Xe+  and  ionic  molecule  Xe~.  (See 


Table  4  for  the  footnotes.) 


N 

Xe  + 

_ n 

Structure  ^ 

R  2> 

I 

Xe2 

3.22 

II 

V 

Ae^ 

linear 

3.383 

III 

X< 

regular  pyramid 

3.646 

IV 

X< 

(Xe^ )  (Xe) 

3'zv  'xz 

3.384 

V 

Xet 

<4 

linear 

3.315 

VI 

Xe5 

(XeJ)ni(Xe)xz 

3.646 

VII 

Xes 

(Xe3>z<Xe2>xz 

3.383 

VIII 

Xe6 

(Xe4> IIIXe2 

3 . 644 

IX 

Xe6 

(Xet)  (Xe-,) 

3'zv  3'XZ 

3.375 

R  3> 

4) 

q 

E 

Dissociation 

Products 

De 

-1.080 

Xe+  +  Xe 

1.08 

-1.277 

Xe£  +  Xe 

0.197 

-1.480 

Xe£  +  2Xe 

0.40 

4.19 

0.014 

-1.356 

Xe^  +  Xe 

0.079 

3.65 

0.184 

-1.346 

Xe^  +  Xe 

0.069 

4.40 

0.001 

-1.568 

Xe£  +  Xe 

0.088 

4.24 

0.005 

- 1.439 

Xe^Xe  +  Xe 

0.083 

4.40 

0 

-1.637 

Xe+Xe  +  Xe 

0.069 

4.20 

0.021 

-1.547 

Xe^Xe2  +  Xe 

0.108 
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Table  6.  Structures  of  the  ionic  clusters  Ar^(ArH)+  and  ionic  molecule  (ArH)+. 


The  distances  (R) 

are 

in  A,  the  energies 

(E)  in  eV 

and  the 

dipole 

moments  (p) 

in  Debyes 

(These  are  taken 

from  Reference  14. 

■  ) 

Ark(ArH)4 

U 

RAr-H 

2) 

R  2) 

Ar-Ar 

M 

E 

Ar 

Detachment 

(ArH) + 

0.569 

2.78 

-4.055 

Ar(ArH)+ 

0.579 

2.81 

3.70 

5.98 

-4.218 

0.163 

Ar2(ArH)+ 

0.591 

2.79 

3.81 

6.32 

-4.395 

0.177 

Ar^ (ArH)+ 

0.601 

2.79 

3.72 

6.23 

-4.586 

0.191 

Ar^(ArH)''" 

0.609 

2.84 

3.57 

4.94 

-4.760 

0.174 

Ar, (ArH) 4 

0.614 

2 . 84 

3. 

17  3.87  3.38 

3.33 

-4.947 

0.094 

^Charge  of  the  H  atom  in  (ArH)  + 

^R,  u  and  R  are  the  distances  from  the  neutral  atom  to  H  and  Ar  of  (ArH)  +  , 

Ar  -  H  Ar-Ar 

respectively.  For  Ar&(ArH)+  there  are  two  different  Ar-H  and  Ar-Ar  distances. 
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Table  7 . 


■jt 

Structure  of  the  Rydberg  excited  Ar2  dimers  and  linear  symmetrical  Ar^ 
clusters  in  metastable  configurations.  The  dissociation  energies  Dg  of 
Ar2  -*•  Ar  +  Ar  and  Ar^  -»  Ar2  +  Ar  atom  detachment  are  in  eV,  and  the 

o 

equilibrium  distances  Rg  are  in  A.  The  distribution  of  the  excitation 
among  the  Ar^  atoms  is  presented  by  q^,  i  -  2  being  the  central  atom. 
(These  are  taken  from  Reference  180.) 


State 

excitation 

energy 


State 

R 


Without  SO 

With  SO 

3P(4s) 

1P(4s) 

1D(4p) 

3P2(4s) 

3Pl(4s) 

1F1(4s) 

11.58 

11.88 

13.09 

11.52 

11.61 

11.87 

X 

u 

A(0u) 

B<°u> 

C<°u> 

2.48 

2.48 

2.46 

2.52 

2.47 

3.60 

0.49 

0.56 

0.79 

0.42 

0.31 

0.055 

2.59 

0.008 

0.16 

0.68 

0.16 


2.58 

0.042 

0.15 

0.70 

0.15 


2.515 

0.232 

0.09 

0.82 

0.09 


3.63 

0.028 

0.31 

0.38 

0.31 
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Figure  Captions 

1.  Ar^  and  Ar^  geometries.  The  Roman  numerals  indicate  the  clusters  (the 
same  as  in  Table  4).  The  numbers  without  signs  are  distant''  in  a, 
and  the  numbers  with  the  +  sign  are  the  atomic  charges.  (This  figure 
is  taken  from  Reference  14.) 

2.  Xe^  and  Xe^  geometries.  For  notations  see  Figure  1.  (This  figure  is 
taken  from  Reference  14.). 

3.  Geometries  of  the  quasistable  XeHCl  clusters  in  the  states  with  electron 

o 

transfer.  The  numbers  stand  for  interatomic  distances  in  A  and  atomic 
charges.  The  energies  (relative  to  the  van  der  Waals  cluster)  of  the 
quasistable  ground-state  isomer  li  and  the  excited-states  isomers  3i  and 
4i  are  4.24  eV,  7.34  eV  and  8.43  eV,  respectively.  (This  figure  is  taken 
from  Reference  39.) 

4.  Energy  spectrum  of  the  Xe^Cl  cage  excitation  from  the  ground  to  excited 
states  with  Xe  -*  Cl  electron  transfer.  The  temperature  is  10  K.  (This 
figure  is  taken  from  Reference  60.) 


